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By splitting a Hamiltonian into two parts, using the solvability of eigenvalue problem of one part of 
the Hamiltonian, proving a useful identity and deducing an expansion formula of power of operator 
binomials, we obtain an explicit and general form of time evolution operator in the representation 
of solvable part of the Hamiltonian. Further we find out an exact solution of Schrodinger equation 
in a general time-independent quantum system, and write down its concrete form when the solvable 
part of this Hamiltonian is taken as the kinetic energy term. Comparing our exact solution with 
the usual perturbation theory makes some features and significance of our solution clear. Moreover, 
through deriving out the improved forms of the zeroth, first, second and third order perturbed 
\^ , solutions including the partial contributions from the higher order even all order approximations, 

. we obtain the improved transition probability. In special, we propose the revised Fermi's golden 

' rule. Then we apply our scheme to obtain the improved forms of perturbed energy and perturbed 

^ state. In addition, we study an easy understanding example to illustrate our scheme and show 

its advantage. All of this implies the physical reasons and evidences why our exact solution and 
O ' perturbative scheme are formally explicit, actually calculable, operationally efficient, conclusively 

I more accurate. Therefore our exact solution and perturbative scheme can be thought of theoretical 

, developments of quantum dynamics. Further applications of our results in quantum theory can be 

expected. 
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One of the most important tasks of physics is to obtain the time evolution law and its form of a physical system, that 
is, so-called dynamical equation and its solution. In a quantum system, the dynamical equation is the Schrodinger 
equation ,1, 2] for a pure state or the von Neumann equation for a mixed state P|. Suppose that, at time t = 0, an 
arbitrary initial state is |5'(0)}, then at a given time t, the state of quantum system will become 

g: \^{t))^e-'"'\^m. (1) 

^ ■ where H is the Hamiltonian of a quantum system and it is assumed independent of time, while e~'^* [h = 1) is 
i~ \ called as the time evolution operator. Above equation describes an arbitrary initial state how to evolve to a final 
^ . state, that is, it is a formal solution of the Schrodinger equation for a pure state. Unfortunately, only finite several 
" quantum systems are exactly solvable at present. This implies that in a general quantum system, the time evolution 
of expanding coefficients of final state in some given representation can not be clearly and explicitly expressed by the 
c-number function. From our point of view, this shortcoming leads to that we can not fully determine and exactly 
calculate those physical quantities dependent on the state (that is expanding coefficients in a given representation) 
at any time. In other words, above formal solution of final state in the operator form is not really useful for the 
many practical purposes. For example, in the cases with environment-system interaction, how does the final state 
decohere, and in the cases with an entangled initial state, how does the entanglement of final state vary, and so on? 

In order to overcome above shortcoming appearing in a general quantum system, quantum perturbation theory 
is well developed and it is successfully applied the calculations of perturbed energy, perturbed state and transition 
probability etc. In the usual perturbation theory, the key idea to research the time evolution of system is to split the 
Hamiltonian of system into two parts, that is 

H = Ho + Hi, (2) 

where the eigenvalue problem of Hq is solvable, and Hi is the rest part of the Hamiltonian except for Hq. In other 
words, this splitting is chosen in such a manner that the solutions of Hq are known as 

/fo|$^)=i?7l*^): (3) 
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where is the eigenvector of Hq and Ej is the corresponding eigenvalue. However, when Hi is not so small compared 
with Hq that we have to consider the higher order approximations or when the partial contributions from higher order 
even all order approximations become relatively important to the studied problems, the usual perturbation theory 
might be difficult to calculate to an enough precision in an effective way, even not feasible practically since the lower 
approximation might break the physical symmetries and/or constraints. In our point of view, the physical reasons 
resulting in above difficulties and obstacles are that the usual perturbation theory does not give a general term form of 
expanding coefficient evolution with time for any order approximation and does not explicitly express it as a c-number 
function. Thus, it is necessary to find the perturbed solution from the low to the high order approximation step by 
step up to some order approximation for a needed precision. These problems motivate us to focus on finding an 
explicit expression of general term for any order approximation, that is an exact solution that is not in the operator 
form, and then reasonably and physically deduce out a scheme of perturbation theory. In fact, as soon as this purpose 
is arrived at, we will obtain an exact solution with all order contributions coming from Hi. This consequentially 
implies, in principle, that this scheme can more effectively and more accurately solve the physical problems than 
using the existing method. Then, we concretely realize and demonstrate this physical idea by providing a scheme of 
perturbation theory. In this process, just like the accustomed way to study a more complicated quantum system, if 
Hi is indeed small enough compared with Hq, its inffuence may be as a perturbation and one can take up to some 
order approximation conveniently at a needed precision. However, the distinct difference between ours and usual one 
is that our perturbed solution will be able to include the partial contributions from the higher order even all order 
approximations based on the reasonable physical considerations so that the physical problems are more accurately 
and effectively calculated, and kept their symmetries and/or constraints as possible since we have used the feature 
of our exact solution with the general term form of time evolution of arbitrary initial state and added systematically 
the contributions from the high order approximations. 

Actually, one of interesting and important tasks in the study of quantum dynamics is to find clearly the the expand- 
ing coefficients (amplitudes) before the basis vectors in a given representation how to evolve. The typical examples 
are the calculations decoherence and disentanglement of quantum systems. Although we split the Hamiltonian into 
two parts in terms of the known method, the problem has not yet been finally answered because that Hq does not 
commute with Hi in general, and the matrix elements of time evolution operator e~'^* do not have the closed factor 
forms about time t in the given representations including Hq representation spanned by {!$'''), 7 = 1, 2, • • • } in terms 
of the known perturbation theory. To solve this kind of problems is also one of our main motivations. 

Just well-known, quantum dynamics and its perturbation theory are correct and have been sufficiently studied. 
Many famous physicists created their nice formulism and obtained some marvelous results. An attempt to improve its 
part content or increase some new content as well as some new methods must be very difficult in their realizations. This 
results in our this possibly lengthy manuscript. However, our endeavor is only at the beginning. It must be pointed 
out we still use the extensively accepted principles and laws of quantum mechanics, and we try to use elementary 
mathematics as possible. All of mathematics knowledge used here does not exceed the university level. 

Undoubtedly, the formalization of physical theory including quantum dynamics often has its highly mathematical 
focus, but this can not cover its real motivations, potential applications and related conclusions in physics. We study 
the reexpression in form should not be only a mathematical game and skill, we expect to find its physical content 
and applications. This is leaded actually by our more explicit and general form of the Schrodinger equation and our 
perturbative scheme. A series of conclusions done here make us assure that we have partially arrived at our purpose. 

In this paper, we start from finding an explicitly closed form about time t of matrix elements of time evolution 
operator e^'-^* in the Hq representation spanned by {!$'>'), 7 = 1, 2, • • • }. Our method is first to derive out an expansion 
formula of power of operator binomials and then apply it to the Taylor's expansion of time evolution operator e^'^*. 
Moreover, by proving and using our identity, we derive out the explicit and general form of representation matrix of 
Q-iHt ^YiQ representation of Hq. Consequently, we obtain an exact solution of the Schrodinger equation in a general 
time-independent quantum system, in particular, its concrete form when the solvable part of Hamiltonian is taken as 
the kinetic energy term. Furthermore, by comparing our solution with the usual perturbation theory, we reveal their 
relation and show what is more as well as what is different in our exact solution. The conclusions clearly indicate 
that our exact solution is consistent with the usual perturbation theory at any order approximation, but also in 
our exact solution we explicitly calculate out the expanding coefficients of unperturbed state in Lippmann-Schwinger 
equation for the time-independent (stationary) perturbation theory and/or fully solves the recurrence equation of 
the expansion coefficients of final state in Hq representation from a view of time-dependent (dynamical) perturbation 
theory. After introducing the perturbative method, in order to provide a new scheme of perturbation theory based on 
our exact solution, we first propose two useful skills. Then, we expressly demonstrate our solution is not only formally 
explicit, but also actually more accurate and effective via generally deriving out the improved forms of the zeroth, 
first, second and third order approximation of perturbed solution including partial contributions from the higher order 
even all order approximations, finding the improved transition probability, specially, the revised Fermi's golden rule, 
and providing a operational scheme to calculate the perturbed energy and perturbed state. Furthermore, by studying 
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a concrete example of two state system, we illustrate clearly that our solution is more efficient and more accurate than 
the usual perturbative method. In short, our exact solution and perturbative scheme are formally explicit, actually 
calculable, operationally efficient, conclusively more accurate (to the needed precision). 

It is worth emphasizing that we obtain an explicit and general expression of time evolution operator that is a 
summation over all order contributions from the rest part of Hamiltonian except for the solvable part, or all order 
approximations of the perturbative part of the Hamiltonian. This makes us be able to build our perturbative scheme 
via so-called "dynamical rearrangement and summation" technology, which is seen in Sec. IVIII In a sentence, our 
perturbative scheme comes from the our exact solution, and then show our exact solution indeed contain useful and 
interesting physical content besides its explicit mathematics form. 

Based on the above all of reasons, our exact solution and perturbative scheme can be thought of theoretical 
developments of quantum dynamics. In our point of view, it is helpful for understanding the dynamical behavior and 
related subjects of quantum systems in theory. Specially, we think that the features and advantages of our solution 
can not be fully revealed only by the perturbative method, because it is an exact solution. We would like to study 
the possible applications to the formulation of quantum dynamics in the near future. 

This paper consists of 11 sections and one appendix. Besides Sec. 13 is an introduction, Sec. IXII is the conclusion 
and discussion, the other 9 sections can be divided into three parts. The first part made of sections from two to four 
is to deduce and prove our exact solution in the general time-independent quantum system; the second part made 
of section five is to compare our exact solution including all of order approximations with the usual perturbation 
theory and demonstrate the features and advantages of our exact solution; the third part made of sections from six 
to ten is to propose our perturbative scheme, obtain its application and illustrate its results. More concretely, they 
are organized as the following: in Sec. [H] we first propose the expansion formula of power of operator binomials; 
in Sec. IIIII we derive out an explicit and general expression of time evolution operator by proving and using our 
identity; in Sec. IIVI we obtain the exact solution of the Schrodinger equation and present a concrete example when 
the solvable part of Hamiltonian is taken as the kinetic energy term; in Sec. 0we compare our solution with the 
usual perturbation theory and prove their consistency, their relations and explain what is more as well as what is 
different in our solution; in Sec. IVII we introduce two skills to relate our exact solution and the perturbative scheme 
in order to include the contributions from the high order approximations; in Sec. IVIII we deduce out the improved 
perturbed solution of dynamics including partial contributions from the higher even all order approximations; in Sec. 
I Villi we obtain the improved transition probability, specially, the revised Fermi's golden rule. In Sec. IIXI we provide 
a scheme to calculate the perturbed energy and the perturbed state; in Sec. ^ we study an example of two state 
system in order to concretely illustrate our solution to be more effective and more accurate than the usual method; 
in Sec. IXII we summarize our conclusions and give some discussions. Finally, we write an appendix where the proof of 
our identity is presented and some expressions are calculated in order to derive out our improved forms of perturbed 
solutions. 



II. EXPANSION FORMULA OF POWER OF OPERATOR BINOMIALS 



In order to obtain the explicitly exact solution of the Schrodinger equation in a general time-independent quantum 
system, let us start with the derivation of expansion formula of power of operator binomials. Without loss of generality, 
we are always able to write the power of operator binomials as two parts 

(yl + B)" = + (4) 

where A and B are two operators and do not commute with each other in general. While the second part /"(A, B) is 
a polynomials including at least first power of B and at most nth power of B in its every term. Thus, a general term 

with Ith power of B has the form (^n'=i A'^^ ^"^'^^^=1 From the symmetry of power of binomials we conclude 
that every ki take the values from to {n — l), but it must keep n — I — — 0- So we have 

(5) 
(6) 



n n- 



1=1 ki,--- ,fc,=0 \«=1 / 
n n — l /I \ 

E E n^'^'^M^ 

1=1 fci,--- ,fci=0 \i=l J 
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where 9{x) is a step function, that is, 9{x) = 1 if a; > 0, and 9{x) = if a; < 0. Obviously based on above definition, 
we easily verify 

f\A,B)^B, f{A,B)=AB + B{A + B). (7) 

They imply that the expression (0) is correct for n = 1, 2. 

Now we use the mathematical induction to proof the expression (|SJ) of /"(v4, B), that is, let us assume that it is 
valid for a given n, and then prove that it is valid too for n + 1. Denoting J-'^^^{A, B) with the form of expression 
© where n is replaced by n + 1, and we extract the part of / = 1 in its finite summation for I 

n n+l {n+l)-l / I \ 

r'+\A, B) = BA''-''' n ^ '=1 

fci=0 i=2 fci,--- ,fci=0 \i=l / 



1) 



(8) 



We extract the terms fci = in the first and second summations, and again replace ki by fci — 1 in the summation 
(the summations for ki from 1 to (n + 1) — / change as one from to n ~ I), the result is 

n n+l {n+l)-l / I \ 

Jf^^+i (A, B) = BA" + Y BA"--^^ + ^ \\\A^^b\ 

fc,=l 1=2 ko.--- .ki=n\i=2 ) 



n+l(n+l)-i (n+l)-i / I \ 
1=2 fci = l k-2.,--- .ki=^ \i=\- / 



(9) 



furthermore 



n-l n+l (n+l)-i / I \ 

fci=0 Z=2 fc2,--- ,fci=0 \i=2 / 



-J^(n+\)-l-Y,\^^k 

n+l n-l (n+l)-( / ; 



e (ii^'''B\A--^-^'^--^'^'e[n-i-YkA. (10) 

;=2 A;i=0 fe2,- - ,fcl=0 \i=l / V i=l / 

Considering the third term in above expression, we change the dummy index {^2, ^3, • • ■ , h} into {fci, ^2, • ■ ■ , ^z-i}, 
rewrite (n + l) — / asn — (/ — 1), and finally replace Z by Z — 1 in the summation (the summation for I from 2 to n + 1 
changes as one from 1 to n), we obtain 

n+l (n+l)-i / I \ / ' \ 

^E E n^''^p("+i)-'-sU'=.0 („+!)_/ 

1=2 k2,--- ,ki=0 \i=2 ) \ i=2 ) 

n+l n-(i-l) n-\ \ / \ 

= ^E E n^'"^U"-('-i)-5:'=J'='(? n-(?-l)-^fc, 

i=2 fci, -- ,fci_i=0 \i=l / \ i=l / 



n n 



-I 



^E E n^''^M""'"^"^^''M^^-^-E^O 

1 = 1 ki,--- ,ki=0 \i=l ) \ i=l / 

BriA,B), (11) 



5 



where we have used the expression © of /"(A, B). Because of the step function 9 {n ~ I ~ X)i=i ^i) ™ t^*^ fourth 

term of the expression Hl()|l . the upper bound of summation for I is abated to n, and the upper bound of summation 
for fc2, fcs, • ■ • ,ki is abated to n — I. Then, merging the fourth term and the second term in ea. (|10|l gives 

n-l n+l n-l (n+l)-l / / \ 

A:i=0 Z=2 fei=0 fc2,--- ,fcj=0 \i=l / 



n 71 — 1 



ay ^'''BA''-^-''' +ay y n^'^'*^ 



n n — l 



Yh 

1=1 



= ^E E [l[A'''B\A^-'-^'-^'''ein 
;=i fei,fe2,- - ,fc,=o \j=i / V 
= Ar{A,B), 

where we have used again the expression © of /""{A, B). 

Substituting pifl and \i2\ into pUf) . immediately leads to the following result 



i-Y^^^ 



i=l 



Note that 



we have the relation 



= BA" + {A + B)f\A, B). 
+ = {A + B){A + BY 

B) = BA"" + {A + B)P[A, B). 



(12) 



(13) 



(14) 



(15) 



Therefore, in terms of eas. p3|l and ()15|l we have finished our proof that /"(A, i?) has the expression ©or © for 
any n. 

III. EXPRESSION OF THE TIME EVOLUTION OPERATOR 



Now we investigate the expression of the time evolution operator by means of above expansion formula of power of 
operator binomials, that is, we write 



(16) 



ri=0 



Tl = 



In above equation, inserting the complete relation |$'')($'''| — 1 before every -ffpS and using the eigen equation 
of iJo ©, it is easy to see that 



r(i/o,i?r) = E E E 



i=l 7i,--- ,71 + 1 fci....,fc,=0 Li=l 

Y.\^^ ki+l<n 



XE. 



."-Ei=i ki-i 

7! + l 



$71 1) ^$71 + 1 I 



Y Y cr{E[j,i]) 

1=1 7l,--- ,7! + l 



J^^7,7.+ l 



|$7l)($7< + l|^ 



(17) 
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where if^'-^'+i 



^$7i [iJj^l^T'+i), E['y,l] is a vector with / + 1 components denoted by 
and we introduce the definition of Cf(i?[7, Z]) (/ > 1) as the following 



E 



fci,-- - ,fc,=0 



.1=1 



E. 



-ELifc.-' 



7! + l 



(18) 



(19) 



In order to derive out an explicit and useful expression of C"(i?[7,Z]), we first change the dummy index ki — > 
n — l — Yl\=i ill tlis eq.®. Note that at the same time in spite of 6{n — / — Yl\=i ^i) changes as 9{ki), but a hiding 
factor 9{ki) becomes 6{n — I — ^i)- Thus, we can rewrite 

1=1 ki,--- ,fci=0 \i=l ) 



(20) 



In fact, this new expression of /"(A, i?) is a result of the symmetry of power of binomials for its every factor. Similarly, 
in terms of the the above method to obtain the definition of C"(£'[7, V\) p9|) . we have its new definition (where fc; is 
replaced by fc) 



n — l\ 71 — 1 

fc=0 I fei,--- ,fc,_i=0 



7-1 



1=1 



(21) 



Because that 9 (^{n — k) ~ I ~ X]'=i abates the upper bound of summation for ki{i = 1, 2, • • • , ^ — 1) from {n — /) 
to (n — fc — 1) — (Z — 1), we obtain the recurrence equation 



A ( (n-fc-l)-(i-l) 



cr(i?[7,z]) = E E 

A:=0 [ fci,--- ,fc,_i=0 



7-1 



,(n-fc-i)-5:^z;fei-(i-i) 

'7! 



1=1 

x0 (^(n - fc - 1) - (/ - 1) - ^ fc,^ I E\^^ 

Y,crs,'-'{E[^,i-i])E';^^^. 



(22) 



fe=0 

In particular, when Z — 1, from the definition of C'"{E{'-f, 1)) it follows that 

n-l 



cr(i?(7, 1)) = E = E 

A:i=0 fci=0 



n-l / j-1 X ''I 
^71 



72 



By means of the summation formula of a geometric series, we find that 



c^r(i^(7,i)) = 



e: 



71 



rpn 
^72 



Eyi -E^72 ^71 -^72 



(23) 



(24) 
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Based on the recurrence equation 122() . we have 



"-2 / 



73 



j._Q \ ^71 ^72 ^71 ^72 
»=1 -^Ti ~ ^■^■^ k=0 ^ -^T- 

^11 ^12 



E 



(^71 ^ ^72) (-^71 ~ ^73) (^71 ~ ^72) (^72 " ^73) 



rpn 

^73 



(-^71 -^73 ) (-^72 -^73 ) 



(25) 



In the above calculations, the last step is important. In fact, in order to obtain the concrete expression of C" {I and 
n are both positive integers), we need our identity 

'+1 f (If < < 

It is proved in Appendix A in detail. The denominators di{E['y, I]) in above identity are defined by 

I 



di{E[^A) = 11(^71 -^^7.+i), (27) 
1=1 

d^{E[^, I]) = \{ {E,, - i?7.) n (^7. -E,,)^ (28) 

j = l fc=i+l 

di+i{Eb,l]) = n (^7. - ^^7,+ J , (29) 
where 2 < i < I. Then, using the recurrence equation H22|) and our identity (|25f) . we obtain 

cr(i=h.,i)^i:(-i)'-'5^, (30, 



1=1 



Here, the mathematical induction shows its power again. If the expression (|30|l is correct for a given n, for example 
n = 1,2, then for n + 1 from the recurrence equation (|22|l it follows that 



n+l-l / I pn-k \ 



. . d,(EhA-l\) ^ 



^71 + 1 



^d.(i?[7,Z-l]) fr^g V ^7, 

^ ' (-irwiz^^:^^ 

^ d,(i?[7,z-i])(£;7.-£^7,+J ■ ^ ^ 

Because that di{E[^, I — l])(i?^. — E^^^^) — di{E[j, I]) {i < I + 1) and our identity H26|) . it is easy to see 



8 



Substitute it into ea. H31() yields 



ly 



(33) 



That is, we have proved that the expression H^Ofl of C"(£^[7, 1]) is vahd for any n. 

It must be emphasized that for the diagonal elements or in the degeneration cases, we need to understand above 



For instance, for C"(i?[7,l]) 



-1 



/ {Ej-^ — Ej^ ) , we have that 



expressions in the sense of limitations. 

the expression lim^e^^^B^^ C^{E[j, 1]) = {Sj^^^ + ©(172 - 1i\)Se-,^e^J nE!^~^, where the step function Q{x) = 1, if 
X > 0, and 8(2:) = 0, if cc < 0. Because of our identity H26|l . the summation to I in ea. p7|l can be extended to 00. 
Thus, the expression of the time evolution operator is changed to a summation according to the order (or power) of 
the Hi as the following 



Ht 



i+i 



lE-y.t 



1=1 7l:--- .7i + l Li=l 

It is clear that it has the closed time evolution factors. 



d,{E[^,l\ 



X{Hr'^^6,,,5,,^,,,. (34) 



IV. SOLUTION OF THE SCHRODINGER EQUATION 



Substituting the our expression of the time evolution operator 134(1 into eq.l^, we obtain immediately our explicit 
form of time evolution of an arbitrary initial state in a general quantum system 



Ar'{t) 



E 

7,7' 



E 



i+i 



71,--- >7l + l 



d,{E[^,l]) 



J|^7,7. + i 

J = l 



'^7l7'^7i + l7' 



(35) 
(36) 



Of course, since the linearity of the evolution operator and the completeness of eigenvector set of Hq, we can simply set 
the initial state as an eigenvector of Hq- Up to now, everything is exact and no any approximation enters. Therefore, it 
is an exact solution in spite of its form is an infinity series. In other words, it exactly includes all order approximations 
of Hi . Only when Hi is taken as a perturbation, it can be cut-off based on the needed precision. From a view of 
formalized theory, it is explicit and general, but is not compact. Moreover, if the convergence is guaranteed, it is 
strict since its general term is known. Usually, to a practical purpose, if only including the finite (often low) order 
approximation of Hi, above expression is cut-off to the finite terms and becomes a perturbed solution. 

Actually, since our solution includes all order contributions from the rest part of Hamiltonian except for the solvable 
part Hq, or all order approximations of the perturbative part Hi of Hamiltonian, it is not very important whether 
Hi is (relatively) large or small when compared with Hq. In principle, for a general quantum system with the normal 
form Hamiltonian, we always can write down 



fr2 

H = — 
2to 



fc2 



V. 



If taking the solvable kinetic energy term as Hq = - — and the potential energy part as Hi 

2m 

applicable to such quantum system. It is clear that 



(37) 



V, our method is 



2m 



|fc) 



2m 



\k) = Ek\k), Ek^ 



2m 



And denoting 



{x\k) 

V'fe(o) 



1 



L3/2 

(fc|*(o)). 



(38) 



(39) 

(40) 
(41) 
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we obtain the final state as 



L3/2 



■J~J ykjkj + i 



EE E '^-(o) 

L3/2 



'5felfc'5fci+lfc' 



(42) 



k 



EE E '^-(o) 

E (~^)' ^ ^ i(fc-n;-Efc.t) 

/n r,3/2 



nV(fc,-fc, + i) 



In the second equal mark, we have used the fact that V = V{x) usually, thus 
^fe.fc.+i = {kj\V{x)\k,+,) 

/oc 
d^Xjd^Xj+i{kj\xj){xj\V{x)\xj+i){xj+i\kj+i) 
-OO 



(43) 



= — / d^x,d^x,+ie-"'^''^Vix,+i)S^ix,-x,+,)e"'^+'-^+'^ 



1 r°° 

= — d=^a;y(a;)e"'('=^-'=^+i) 

L J-oo 

= V{kj - kj+i). 



(44) 

It is clear that V(fe) is the Fourier transformation of V{x). The ea. l42() or (|43() is a concrete form of our solution 
(|35|l of Schrodinger equation when the solvable part of Hamiltonian is taken as the kinetic energy term. Therefore, 
in principle, if the Fourier transformation of V{x) can be found, the evolution of the arbitrarily initial state with 
time can be obtained. It must be emphasized that it is often that the solvable part of Hamiltonian will be not only 

for the practical purposes. Moreover, since the general term is obtained, we can selectively include the 



T 



2m 



partial contributions from some high order even all order approximations. 
Furthermore, we can derive out the propagator 



(a;|e~'^*|a;') = {x\k){k\e-'"'\k'){k'\x') 

kM' 

oc 

= Eia^'^^-^'^-'^^+EE E 



l — l k,k' fci ,■ ■ ■ .k 



■J~J ykjkj + i 



5kik5ki + ik' 



{-ly 



^^d,{E[k,l]) L3/2 



^ifc-a; — ifc -x 



L3^ 



(45) 



It is different from the known propagator in the expressive form, but they should be equivalent to the physical results. 
One of its distinguished features is this matrix element to be fully c-number function made of the matrix elements 
V'^^'^^ and unperturbed energy levels E^^ . 



V. COMPARING WITH THE USUAL PERTURBATION THEORY 



In this section, we will compare the usual perturbation theory with our solution, reveal their consistency and 
relation, and point out what is more in our solution and what is different among them. Furthermore, we expect to 
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reveal the features, significance and possible applications of our solution in theory. We will respectively investigate 
and discuss the cases comparing with the time- independent, time-dependent perturbation theories as well as the 
non-perturbed solution. 

A. Comparing with the time-independent perturbation theory 

The usual time-independent (stationary) perturbation theory is mainly to study the stationary wave equation in 
order to obtain the perturbative energies and perturbative states. However, our solution focus on the development of 
quantum states, which is a solution of the Schrodinger dynamical equation. Although their main purposes are different 
at their start point, but they are consistent. Moreover, our exact solution is also able to apply to the calculation of 
perturbed energy and perturbed state, which will be seen clearly in Sec. IIXI 

In order to verify that the time-independent perturbation theory is consistent with our solution, we suppose that 
in the initial state, the system is in the eigen state of the total Hamiltonian, that is 



(46) 



It is clear that 



7,7' I ' = 1 J 

To calculate the if th time derivative of this equation and then set t = 0, we obtain 



(47) 



7,7' L 1 = 1 



where we have used the fact that 



E <v+E^7''wr«7'i*^), 

7,7' I 



1 = 1 



I'^eAo)) - E«7l*^), «7 = {^'\^eAo)), 



(48) 
(49) 

(50) 



d^A]-' it) 



dt^ 



t=o 



(51) 



It is easy to obtain that 

d^A7^'(t) 



dt^ 



= (-i)" E 



t=0 



7l,--- •ll + l Li=l 



l+l pK 

E(-ir^^^ 



'^7l7'^7i + l7'i 



71 ,11 + 1 

where we have used our identity (|26|l . This means that iil > K 



-Q^7,7, + l 



'^7i7'^7l + l7' ■ 



(52) 
(53) 



di^ 



{~i)''Br'{K>l)^0. 



t=o 



In special 



(54) 



it) 

dt 



-isr (1) = -mr sn. 



(55) 
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A^Ar it) 



di2 



(56) 



In other words, the summation to I in the right side of ea. H48f) is cut-off to K. Therefore, by left multiplying {^''\ to 
ea. (|48|l we obtain 



K 



7'- 



7' 1=1 

When K = 1, above equation backs to the start point of time-independent perturbation theory 



(57) 



(58) 



It implies that it is consistent with our solution. In order to certify that our solution is compatible with the time- 
independent perturbation theory, we have to consider the cases when K > 2. Note that Ei^a^ = Et [E^^^ 
substitute ea. H57|l twice, then move the second term to the right side, we have 



K 



K-l 



E^-'Era, = <a, -I- J] J] i?7^'(i^)ay - ^ J] i?,''^'(i^ - 1) J] (i?y<5y,. + H^'^'') a 



7' 1 = 1 

K-l 



7' 1=1 



l" 



<«7 + E E [bv\k) - EyBr\K - 1) 



7' 1 = 1 
K 



K-1 



Since when K > 2 



and again / > 2 



E n^i^'"^ ^77.^7,+.7'«7' -E E Br\K-i)Y,Hr"c 

7' \i=l / 7' '=1 7" 



Cf (i?[7, 1]) - i?7,Cf -i(i?[7, 1]) = <- 



C^{E[^,l]) ~ E,,^,C^-\E[^,l]) = Cl-\E[^,l- 1]), 
which has been proved in appendix A. Obviously, when K — 2, ea. (|59|l becomes 

E^Etcl^ = E'^Qj -|- E^ ^ ^ H^"^ a.yi . 



7' 



It can back to ea. H58(l . Similarly, for X > 3, so do they. In fact, from eas. (|6()l61|l . ea. (|59|l becomes 

E^-^Era^ = ^^a^ + E ^7^'«7'' 



7' 



where we have used the fact 

K~l 



E E [^r'iK) - E,.B-p\K - 1)] ^7' + E ( n ^1'"'"' 



'^77l'^7l + l7"^7' 



7' 1=2 

K-1 



7' \J = 1 



EE^'''(^-i)E<'V = o. 



a^], we 



T "7 7 1 



(59) 

(60) 
(61) 
(62) 

(63) 



(64) 
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Its proof is not difficult because we can derive out 

K-l 



J\ — 1 



{K)-E,,Br (K-l) 



f' 1=2 

7' 1=2 7i,- - ,7!+l V = ^ I 

= E E E Cl-\EYi, ^ - 1]) ( ri ^l'"'" ) '577.^7, + .7'«7' 
7' i=2 7i,--- ,7i + l Vj^I / 

= EE E cr\E\iA)\iiiir'Aiir'''^...a,. 

7' i = l 7l>---.7i + l y=i / 

-E E^^''''(^- (65) 

7^,7' l=\ 

The first equality has used the definitions of B'J'^ (K), the second equality has used eq-JHU, the third equality sets 

I — 1 ^ I and sums over 7;+2, the forth equality sets 7' 7", inserts ^7,^17' and uses the definition of B^'^ (K) 

again. Substituting it into the left side of eq.ljHU and using B'^j^_^{K — 1) expression, we can finish the proof of 
ea. (lfi4ll . Actually, above proof further verify the correctness of our solution from the usual perturbation theory. 

In our point of view, only if the expression of the high order approximation has been obtained, we can consider its 
contribution in the time-independent perturbation theory since its method is to find exactly a given order approxi- 
mation. However, this task needs to solve a simultaneous equation system, it will be heavy for the enough high order 
approximation. In Sec. IVIII we will give a method to find the improved forms of perturbed energy and perturbed 
state including partially high order even all order approximations. 

Now, let us see what is more in our solution than the nonperturbative method (sometime be called time-independent 
theory). Actually, since H is not explicitly time-dependent, its eigenvectors can be given formally by so-called 
Lippmann-Schwinger equations as follows f3| : 



|*^(±)) = + -———H,\%{±)) (66) 



1 

E-^ - -Ho ±17?' 

= \^-^)+Gl{±)H^\^l{±)) (67) 

= + W-^^il'f^) (68) 

bj^ — H ±iri 

= 1$'') +G^(±)ffi|$'^}, (69) 

where the complete Green's function G'*'(±) — l/{Ey — H ±\rf) and the unperturbed Green's function Gq(±) = 
1/{E'^ — Hq ± ir;) satisfy the Dyson's equation 

G^(±) = G2(±) + G2(±)i/iG^(±) (70) 

00 

= ^(G2(±)i/i)'G2(±). (71) 

(=0 

Here, the subscript "S" means the stationary solution: 

H\m±))^E,\^li±)). (72) 

More strictly, we should use the "in" and "out" states to express it 1^. In historical literature, this solution is known 
as so-called non-perturbative one. It has played an important role in the formal scatter theory. 
Back to our attempt, for a given initial state |4'(0)), we have 

l*(0)) = 5](vl'i'(±)|*(0))|<(±)). (73) 

7' 
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Acting the time evolution operator on it, we immediately obtain 
l*W> = E(*s'(±)l*(0))e-'^^'*<(±); 



oo ^ oo 

^ [e-'^-%y($7Xo)) 



$7 



00 ^ 

^ (ci>v|(i/,G2;'(T)) |*(o))($^i (G2;'(±)i?, 



(j)7 



;,;'=o 



(74) 



comparing this result with our solution H35|l . we can say that our solution has finished the calculations of explicit 
form of the expanding coefficients (matrix elements) in the second term of the above equation using our method and 
rearrange the resulting terms according with the power of elements of representation matrix H^-l'^^^^ . It implies that 
our solution can have more and more explicit physical content and significance. Therefore, we think that our solution 
is a new development of the stationary perturbation theory. In spite of its physical results which consistent with the 
Feynman's diagram expansion of Dyson's equation, its form is new, explicit, and convenient to calculate the time 
evolution of states with time. Specially, it will be seen that our perturbative scheme based on our exact solution has 
higher efficiency and higher precision in the calculation from Sees. IVIIII Hxl and 1x1 



B. Comparing with the time-dependent perturbation theory 

The aim of our solution is similar to the time-dependent perturbation theory. But their methods are different. The 
usual time-dependent perturbation theory 2] considers a quantum state initially in the eigenvector of the Hamiltonian 
Hq, and a subsequent evolution of system caused by the application of an explicitly time-dependent potential V = 
Xv{t), where A is enough small to indicate y as a perturbation, that is. 



dt 

with the initial state 



i^\^it)) = {Ho + Xvitmit)) (75) 



l*(0)) = I*"). (76) 



It is often used to study the cases that the system returns to an eigenvector |<I>/) of the Hamiltonian Hq when the 
action of the perturbing potential becomes negligible. 

In order to compare it with our solution, let us first recall the time-dependent perturbation method. Note that the 
state \'^{t)) at time t can be expanded by the complete orthogonal system of the eigenvectors 1$''') of the Hamiltonian 
Hq, that is 

\^{t)) = Y.c,{m^) (77) 

7 

with 

c,it) = (<i>^|*(t)), (78) 

whereas the evolution equation (|75|l will become 

i^c^(t) = E^c^{t) + Xj2y''^\t)cy{t), (79) 
7' 

^ . Now setting 

c^{t) = b^{t)e-'^^* , (80) 



where v"''^' {t) = ($^|w(i) 



and inserting this into (|79|l give 



7' 



Then, making a series expansion of b-y{t) according to the power of A: 

6,(t) = 5]A'6«(t), 

and setting equal the coefficients of A' on the both sides of the equation H81(l . ones find: 

i|6«W = ^e'(^-^.')*.--''W&^'r^Ht), (/^O). 
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(81) 
(82) 

(83) 
(84) 



7' 



From the initial condition (|76|l . it follows that c^(0) = ^7(0) = S^a — bj (0) and bj (0) — if I > I. By integration 
over the variable t this will yield 

6f (t) = S^a, (85) 

&(i)(i) = -i /'diie'(^--^°)*iw^"(ii), (86) 



&f(<) = -i rdt2Ve'(^-'-^-)*^t;^^Ht2)&(^,)(i2) 
•^0 72 

= / di2 /" 'dtie'(-^^-^^2)*2eKS-,2-£;»)*i„772(^2)v^^"(ti). 



72 



Jo 



(87) 



If we take v independent of time, the above method in the time-dependent perturbation theory is still valid. Thus 

1 



— Ea 



1 



72 



{E^ — Ea){Ej^ — Ea 



E 



(-E72 - Ea){Ej - Ey^ ) 



('e~'^^2* _ e"'-®^*) e'-^Tt^772 ^72Q 



E 

72 

+ 



1 



1 



(£'7 — Ej2){E^ — Ea) {Ej — Ey2){Ey2 — Ea 
1 



-e^'^^2*e'"^^* 



{E^ — Ea){E^2 ~ ^a) 



g-iBotgiB^t 



(89) 



This means that 



7172 z— 1 



d^iE[J^]) 



^ 77l^72Ct' 



(90) 



.(2) 



1 

(t) = Y^f-_iy-l p-^^-ii^S S ,,7i72„.7273 



71,72,73 i=l 

Now, we use the mathematical induction to prove 

4'Ho= E E(-ir" 

7l,--- ,7! + l i=l 



d.(i?[7,Z]) 



(91) 



(92) 
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Obviously, we have seen that it is vahd for I = 1,2. Suppose it is also valid for a given n, thus form eas. (|84|l and 
6^'^(0) = 0, ; > 1 it follows that 



b^;+'\t) = -i Tdr Ve'(^''-^-)^i;'3^6^")(T). 

Substitute eas. l|5n|) and we obtain 

b'-;+'\t) = -i /'drVe'^^^t;'''^4")(T) 

^ 7 7l,--- :7re+l 2=1 \.7 = 1 / 



(93) 



77l"7»i+ia 



In the last equality, we have inserted ^7„+2/3j summed over 7 and integral over r. In terms of the definition of 

di{E[j, n], we know that di{E[j, n]){E^. — E^^,^^) = di{E['^, n + 1]). Then based on our identity we have 

n+l 

Vf-1)*"^ = = -('-i)("+2)-i ! f95) 



Thus ea. (|94|l becomes 



.^^ d,(£;[7,n+l]) ^ ^ d„+2(£;[7,ri + l]) 



n+2 



E E(-l)'" ^^(j;[^,„'+l]) I n-^-^-"^ I ^^--"^7.../^^7...^e'^^^ (96) 



7l,--- ,7n+2 i=l \j = l 



Set the index taking turns, that is, 71 ^ 7n+2 and 7^ + 1 — >7i,(ri + l) > i > 1). Again from the definition of (ii(i? [7, rt], 
we can verify easily under above index taking turns, 



"+2 -iE^.t "+2 „-i-B,.t 



'd.(i?[7,n+l]) 'd,(£;[7,n+l])- ^^^^ 

Therefore we obtain the conclusion 

n+2 ( / \ 

This implies that we finish the proof of expression (|92|l . Substitute it into ea. (|77|l . we immediately obtain the same 
expression as the our solution JSHJ. 

From the statement above, it can be seen that our solution actually finish the task to solve the recurrence equation 
of bi^{t) by using our new method. Although the recurrence equation of bly{t) can be solved by integral in principle, 
it is indeed not easy if one does not to know our identity H26|l and relevant relations. At our knowledge, the general 
term form has not ever been obtained up to now. This is an obvious difficult to attempt including high order 
approximations. However, our solution has done it, and so it will make some calculations of perturbation theory more 
convenient and more accurate. 

However, we pay for the price to require that H is not explicitly dependent on time. Obviously our solution is 
written as a symmetric and whole form. Moreover, our method gives up, at least in the sense of formalization, the 
requirement that is a small enough compared wit Hq, and our proof is more general and more strict. In addition, 
we do not need to consider the cases that the perturbing potential is switched at the initial and final time, but the 
perturbing potential is not explicitly dependent of time. Since we obtain the general term, our solution must have 
more applications, it is more efficient and more accurate for the practical applications, because we can selectively 
include partial contributions from high order even all order approximation. This can be called as the improved forms 
of perturbed solutions, which will be given in Sec. IVIII 
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VI. TWO SKILLS IN THE IMPROVED SCHEME OF PERTURBATION THEORY 



For the practical applications, we indeed can take the finite order approximations, that is, a finite I is given 
in our solution (|35|l . Moreover, in order to show the advantages of our exact solution and provide the improved 
scheme of perturbation theory, we will derive out the improved forms of perturbative solution including the partial 
contributions from the higher order approximations. Because all of steps are well-regulated and only technology is to 
find the limitation of primary functions. In other words, the calculation of our exact solution and perturbation theory 
is operational in the practical applications, and further it can advance the efficiency and precision in the calculations. 
Frankly speaking, before we know our exact solution, we are puzzled by too many irregular terms and very trouble 
dependence on previous calculation steps. Moreover, we are often anxious about the result precision in such some 
calculations because those terms proportional to t°e~''^^i* (a = 1, 2, • • • ) in the high order approximation might not 
be ignorable with time increasing. Considering the contributions these terms can obviously improve the precision. 
However, the known perturbation theory does not give the general term, considering this task to add reasonably the 
high order approximations is actually impossible. 

Since our exact solution has given the explicit form of any order approximation, that is a general term of an 
arbitrary order perturbed solution, and their forms are simply the summations of a series. Just enlightened by this 
general term of arbitrary order perturbed solution, we use two skills to develop the perturbation theory, which are 
respectively expressed in the following two subsections. 



A. Redivision of Hamiltonian 



The first skill is to decompose the matrix elements of Hi in the representation of Hq into diagonal part and 
nondiagonal part: 



(99) 



so that the concrete expression of a given order approximation can be easily calculated. Note that h^^ has been chosen 
as its diagonal elements and then g'l'^^^^ has been set as its nondiagonal elements: 



i/l — i/l \^ "lj-f] + l)- 

As examples, for the first order approximation, it is easy to calculate out that 



(100) 



-iE-,.t 



{hi (5^j^2 ) (5^^j 5-y'^2 — . e ' ^ 



1! 



1^5 I 

^77 1 



(101) 



AT [9) = E 



71:72 



9i (-^ ^7172)^771^7 



(1 - 5^^,) 



(102) 



Note that here and after we use the symbol AJ'^ denoting the contribution from ith order approximation, which 
is defined by H36() . while its argument indicates the product form of matrix elements of Hi in Hq representation. 
However, for the second order approximation, since 



JJ^ Hi ^ — {hj ) Sj-^j2^i2i3 ~^ 9i '''^^7172 + 3i ^1 ^7273 ~^ 9i 9i '^^ ■ 
we need to calculate the mixed product of diagonal and nondiagonal elements of Hi. Obviously, we have 



(103) 



Ar'{hh)= E 



71.72,73 



E(-i)'-^ 



.4=1 



{hi 57172^1 ^7273)^771^7 



2! 



(104) 
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71.72,73 

e 



i=l 



+ 



^1 ffl '^7l72'^77l'^7'73 



(105) 



71,72,73 Li=l 



rfi(i^[7,2]) 



^7 ^7 7 ^^^^^J^^^J^ 



AV\99) = E 

71,7 

= E 



57" K , 



(106) 



71,72,73 



d^iE[J,2]) 



7172 7273 1 r 
yi J/i "77i"7'73 



Q—iE^t 



Q—iE^^ t 



i^Ej — -^71 ) (-^7 — E^f ) (-^7 — -^71 ) (-^71 — En^' ) 



+ 



g-iEyt 



{Ej E^')(^Ej^ Ery') 



(107) 



In the usual tinio-dcpendent perturbation theory, the zeroth order approximation of time evolution of quantum state 
keeps its original form 



|*(o)(t)> 



-iE^ti 



$7) 



(108) 



where wc have set the initial state as 1$'') for simplicity. By using our solution, we easily calculate out the contributions 
of all of order approximations from the product of completely diagonal elements h to this zeroth order approximation 



E 



7i,'" ,7i+i 



7+1 



E(-ir^ 



iE^.t 



.i=l 



di{E[j,l]) 



n '''l'^7j7j + i ) '^77i^7'7!+i 



(109) 



Therefore, we can add the contributions of all of order approximation parts from the product of completely diagonal 
elements h to this zeroth order approximation to obtain 

*'(°)(t)^ =e-'K+''7)«|$7). (110) 

Similarly, by calculation, we can deduce out that up to the second approximation, the perturbed solution has the 
following form 



|*'(t)) = ^ J e-K^-+''7)M^y + 

7,7' [_ 

+E 



^-i{E^+hDt _ g-i(Ey+h7 ) 

[E^ + hi) - (e^. + hX) 

Q-i{E-,+hj)t 



[{E^ + hi) - {E^, + hi)] [{E^ + hi) - (e^, + hX) 

^-■,{E-,^+h1^)t 



[{E^ + hi) - {E^, + hi)] [{E^, + hH) - (Ey + hX) 

i(^E^,+hX)t 



+■ 



{E^ + hi) - (iJy + hX)] [{E^, +hl)- (Ey + hX 
(i>T'|*(0))l \^^) + 0{Hl). 



Ill 111 

9i gi 



(111) 
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However, for the higher order approximation, the corresponding calculation is heavy. In fact, it is unnecessary to 
calculate the contributions from those terms with the diagonal elements of Hi since introducing the following skill. 
This is a reason why we omit the relevant calculation details. Here we mention it only for verifying the correctness of 
our exact solution in this way. 

The results (|110ll and (|lll|l are not surprised because of the fact that the Hamiltonian is re-divisible. Actually, we 
can furthermore use a trick of redivision of the Hamiltonian so that the new Hq contains the diagonal part of Hi, 
that is 

H'o = Ho+Y.hjmi^'Tl (112) 

7 

H[ = Hi-j2h'im{'^^\^J29T''m{m im 

7 7>7' 

In other words, without loss of generality, we always can choose that H[ has only the nondiagonal elements in the Hq 
(or Hq) representation and 

- (^^ + /i?)l$T) = (114) 

In spite of our skill is so simple, it seems not be sufficiently transpired and understood from the fact that the 
recent some textbooks of quantum mechanics still remain the diagonal contribution from the unperturbed part in the 
expression of the second order perturbed state. 

It must be emphasized that the skill to redivide Hamiltonian leads to the fact that the new perturbed solution can 
be obtained by the replacement 

E^^ Ej^ + hi' (115) 

in the usual perturbed solution and its conclusions. This implies that there are two equivalent ways to obtain the 
same perturbed solution and its conclusions. One of them is to redefine the energy level E^^ as E'^. , think E'^. to 
be explicitly independent on the perturbed parameter from a redefined view, and then use the method in the usual 
perturbation theory to obtain the result from the redivided H'l . The other way is directly deduce out the perturbed 
solution from the original Hamiltonian by using the standard procedure, but a rearrangement and summation are 
carried out just like above done by us. From mathematical view, this is because the perturbed parameter is only 
formal multiplier and it is introduced after redefining E'^^ . The first way or technology will be again applied to our 
scheme to obtain improved forms of perturbed energy and perturbed state in Sec. IIXI 

For simplicity, in the following, we omit the ' in Hq, Hi as well as E^, and always let Hi have only its nondiagonal 
part unless particular claiming. 



B. Contraction and anti-contraction of nondiagonal element product 

In this subsection, we present the second important skill enlightened by our exact solution, this is, a method to 
calculate the contributions from the contractions and anti-contractions of nondiagonal element product in a given 
order approximation. This method is also the most important technology in our scheme of perturbation theory. 

Let us start with the second order approximation. Since we have taken only with the nondiagonal part 

^7j7j+i^ the contribution from the second order approximation is only A^'^ (gg) in ea. HlU7|) . However, we find that 

the limitation in the expression of A'^'^ (gg) has not been completely found out because we have not excluded the case 
Ej = E^i (or 7 = 7'). This problem can be fixed by introducing a decomposition 

gr^gr' = gr^gr^S,,,, + gr^gr^v,,,,, (116) 

where 777173 = 1 — (^-yj^g. Thus, the contribution from the second order approximation is made of two terms, one 
so-called contraction term with the 5 function and another so-called anti-contraction term with the rj function. It 
must be emphasized that we only consider the non-degenerate case here and after for simplification. While when the 
degeneration happens, two indexes with the same main energy level number will not have the anti-contraction. 

In terms of above skill, wc find that the contribution from the second order approximation is made of the corre- 
sponding contraction- and anti-contraction- terms 



at' {gg) = at' (gg; c) + at' (gg; n) , 



(117) 
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where 



71 >7 

= E 



E(-i)^"^ 



71 



+ 



7172 7273 r r r 
i/l yi "7i73"77i"7 n 



(118) 



'1,7 

= E 



71,72,73 Lj=l 



7172 7273 r r 

ai al '/7i73"77i"7 73 



+ 



— — £''y') (-^7 — -^7i)(-^7i — -^7') 

■57"^57^^%7' 



(119) 



This method can be extended to the higher order approximation by introducing a concept of 5-product decomposition. 
For a sequential product of nondiagonal elements g with the form HfeLi fi'i*''"'^^ ("^ ^ 2), we define its (m — l)th 
decomposition by 



n5r''"^= n^r^'M Wi+ lis? 



7fc7fe+l 



(120) 



When we calculate the contributions from the nth order approximation, we first will carry out n — 1 the first decom- 
positions, that is 



fe=i 



.k=l 



7)c7fe+2 



+ 7*1+2) 



Obviously, from the fact that Hi is usually taken as Hermit one, it follows that 



_737j+i _7j+i7j+2 r 

9l 9l °l3lj + 2 



|„737j + l I 

\9i 



-"7373 + 2 • 



(121) 



(122) 



When considering the contributions a given order approximation, the summation over one of two subscripts will lead 
in the contraction of 5-production. More generally, for the contraction of even number g-production 



■J~j'^73 7j+i (5^^^^_^^ j 6^^^6^^^^Y — \9V^^ 
Vi=i fc=i / 

and for the contraction of odd number (^-production, 

m—l 



(123) 



.k=l 



n^"'"^ n^7.7.+2 



'^7l7'^7m+l7' 



fc=l 



9V 



W <^7fc7fc+2 ) '^7l7'^7m+l7'5l'*^ ) 



(124) 



.fc=i 



where ^7i7^7„+i is a factor appearing in the expression of our solution. 

Then, we consider, in turn, all possible the second decompositions, the third decompositions, and up to the (n — l)th 
decompositions. It must be emphasized that after calculating the contributions from terms of lower decompositions, 
some of terms in the higher decompositions may be trivial because there are some symmetric and complementary 
symmetric indexes in the corresponding result, that is, the product of this result and the given (J^^^^, or rj^kiy zero. 
In other words, such some higher decompositions do not need to be considered. As an example, let us analysis the 
contribution from the third order approximation. It is clear that the first decomposition of a sequential production 
of three nondiagonal elements becomes 



„7l72 _72 73 „7374 _ „7l72„7273 „7374 r r 

9i 9i 9i — 9i 9i 9i "7173 "7274 



I 7172 7273 7374 c 
'91 9l 9l "7l73'/7274 
I 7172 7273 7374 c , 7172 7273 7374 

'9l 9l 9l 'hlj3"j2-f4 ' 9l 9l 9l '/7l73'/72l 



(125) 
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Thus, the related contribution is just divided as 4 terms 

(999) = (999; cc) + {999-, cn) + AP' [ggg- nc) + Ap' [ggg, nn) . 



(126) 



In fact, by calculating we know that the second decomposition of the former three terms do not need to be considered, 
only the second decomposition of the last term is nontrivial. This means that 



(555; nn) = AJ^ {ggg; nn, c) + Aj' {ggg; nn, n). 



(127) 



where we have added 5^^^^ in the definition of A'P {ggg;nn,c), and rj-y^^^ in the definition of A'P {ggg;nn,n). Ob- 
viously, in the practical process, this feature largely simplifies the calculations. It is easy to see that the number of 
all of terms with contractions and anti-contractions is 5. For convenience and clearness, we call the contributions 
from the different terms in the decomposition of ^-product as the contractions and anti-contractions of (7-product. Of 
course, the contraction and anti-contraction refer to the meaning after summation(s) over the subscript(s) in general. 
Moreover, here and after, we first drop the argument gg ■ ■ ■ g m. the ith order approximation Ai since its meaning has 
been indicated by i after the Hamiltonian is redivided. For example, the explicit expressions of all contraction- and 
anti-contraction terms in the third order approximation A^ can be calculated as the following: 



2e-'^^ 



di{E[^,i\ 



W9l 



I3I3+1 



-iE-,t 



—3 + {-\t) 2 

{E^ — -E-y') (-^7 ~ -^7') 



{E^ - Ey) 



77 

9r 



9V 



(128) 



AT (cn) 



E 



g— iS^i t 



7ir-- ,11+1 Li=l 

= E 



di{E[^,i\) 



3 

n 



9i 



'^7173^7274^717*^71+17' 



^—iE^t 



^ — iE~ft 



{E^ — -^71) (-^7 — -^7') (-^7 — -^71) (-^7 — ^7') 



-I- 



^—'lEy^t 



(£//y — -^71) (-^71 — -^7') (-^7 — -^7')^ (-^71 — -^7') 



^7 -^71/ V-^71 -^7 



(£//y -^71) (-^7 -^7') 



Ifl7"^l'ff7"'r/7i7' 



(129) 



^^'(r^c) = E 
71. •• 

= E 



E(-ir^ 



lE^.t 



.1=1 



rfi(£;[7,3]) 



^?7l73^7274'^7l7^7i+l7' 



71 



\E-,t 



{E^ — -^71) (-^7 — -^7') (-^7 — -^71) (^71 — -^7') 



(iz-y — -£'71) (-^7 — -^7') (-^7 — ^71) (^71 — -^7') 



^7 -^71/' V-^7 -^7 



(^//y Eryt ) {Ery^ £/,y/ ) 



7 ^71/ V^7l -^7 
2 

^771 ' 



^7" 



^717 



(130) 
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AT {nn, c) 



E 

7172 



c/i(£[7,3]) 



^717*^747' '77173 '?7274 '^77' 



+ ■ 



{E^ - E^^f [E^^ - E^^ ) {Ej - -£72 )^ {E^i - E^^ ) 



{Ej - E^^f [E^-^ - E^^) {E^ - E^^f {Ej^ - E^^) 



-"1 -^71/ \ 71 72 



(-^7 ~ -^71) (-^71 ~ -^72) 



'7 V-^71 ^11) 

,111 „lll2 „l2l' , 



111 1112 121 c 

9i 9i 9i Vii2"- 



111 



(131) 



Aj'''{nn,n) = ^ 



E(-ir^ 



7l.--- ,74 Li=l 

+ E 



7172 



rfi(£;[7,3]) 



„7j73 + l 



(-£^7 ~ -^71) (-^7 ~ -^72) {^1 - ^i') 





e 






{E, 


- ) {E^^ 


- -E'72) (-£^71 


-Ey) 




e" 






{E^ 


- E^^) {Ej^ 




-Ey) 




e" 


-lE^it 




{E, 


- Ej>) {E^^ 


- Ey) (-E'72 


-Ey) 



^71 7^747' ^71 73 ^7274 ^77' 



„771„7172 727 „ „ ,„ , 

yi yi yi Vii2yiii'Vii' 
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In above calculations, the used technologies are mainly to find the limitation, dummy index changing and summation, 
as well as the replacement gj*^^ Vnij = qT'^^ since g^'^^ has been nondiagonal. 

It must be emphasized that, in our notation, A]'^ represents the contributions from the ith order approximation. 
The other independent variables are divided into i — 1 groups and are arranged sequentially corresponding to the 
order of ^(-product decomposition. That is, the first variable group represents the first decompositions, the second 
variable group represents the second decompositions, and so on. Every variable group is a bit-string made of three 
possible element c, n, k and its length is equal to the number of the related order of g-product decomposition, that is, 
for the jth decompositions in the ith order approximation its length is i — j. In each variable group, c corresponds 
to a (5 function, n corresponds to a ry function and k corresponds to 1 (non-decomposition). Their sequence in the 
bit-string corresponds to the sequence of contraction and/or anti-contraction index string. Prom the above analysis 
and statement, the index string of the jth decompositions in the i order approximation is: 



n (7fe,7fc+i+j)' 



(133) 



fe=i 



For example, the first variable group is cccn, which refers to the first decomposition in five order approximation 
and the contribution term to include the factor (57173 (^7274 (^7375 ^774 75 the definition of Ac,{cccn). Similarly, cncc 
means to insert the factor 57i73??7274<57375<57475 in the definition of ^5(cncc). When there are non trivial second 
contractions, for instance, two variable group {ccnn,kkc) represents that the definition of A5{ccnn,kkc) has the 
factor ( (^71 73 <572 74^773 75*774 715) '57376- Since there are fully trivial contraction (the bit-string is made of only k), we omit 
their related variable group for simplicity. 

Furthermore, we pack up all the contraction- and non-contraction terms in the following way so that we can obtain 
conveniently the improved forms of perturbed solution of dynamics inchiding the partial contributions from the high 
order approximation. We first decompose A'^ , which is a summation of all above terms, into the three parts according 



to e 



(-ii)e- 



and (-it) e 



/2: 



-AT (te) 



(134) 



22 



Secondly, we decompose its every term into three parts according to e '^T^*,e ti* e '^71*) and e 



-iE,t. 



Aj-^'ie) = ^^T'(e-'^^*) + A^^'(e-'^^i*)+^^^'(e-'^v*), 
Aj'^'ite) = ^JT'(te-'^^*) + Ar'(te-'^^i*)+^^^'(ie-'-^v*), 



(135) 
(136) 

AJ^' (t^e) = AJ''' (t^e-'^^*) + AJ''' (i^e"'^^!*) + AJ''' (i^e-'^v*). (137) 
Finally, we again decompose every term in above equations into the diagonal and non-diagonal parts about 7 and 7': 



^77(g-iB^,t) ^ (e-'^7^*;D)+A^^ (e-'^-.*;N), 
A^'''(tV^^^*) = yl^^'(i2e-'^-i*;D)+yl^^'(i2e-'^7i*;N), 



where E^. takes , iJ-yj and E^i . 

According above way, it is easy to obtain 



AT (e- 



;D) = - E 



+ - 



1 



{E^ - EjJ {Ej - E^J 



{Ej - EjJ {E^ - E^^) 

^771 „7l72„727J 

91 9i 9i O77') 



AT (e-'^-*; 



N) 



7 ^71 

E-yt 



1 



+- 



(S-y -2'7l) (-^7 -^7') 



„77l „7l7„77 

Si yi yi 



(£/,y -^71 ) (-^7 



{Ej - Ej^) {E^ - Ej^) {E^ - Ey) ' 



AV (e 



^77 (e-iis,it 



'^71*; D) 
N) 



71.72 



{E^ - E^^) (E^i - E^^) 



iE-y^ t 



^771^7172^727 j^^^^,r]^^, 



„77l „7l72 „727X 
^ g-iE^, t 5l ffl ffl 0^7' 

72 

7ii72 

E 

71, 7s 



^77 (e-'^72*;D) 



(-£^7 ~ -^71) (-^71 ~ -^72) (-2^71 ~ ^I'Y 

„77l „7l72„727A , 

i/i yi ill "77 



(-^7 ~ ^72) (-^71 ~ -^72) 



^77 gj gj '?772 ^77' 



71)72 



A^^'(e-'^7'*;N) 



= 0, 



(-E'7 ~ -^72) (^'71 ~ -^72) (-E'72 - Ey) 



1 



+- 



1 



(£//Y — -^7') (-^71 — -E/'-y' 



„7 7l„7i7 „77 

yi yi yi 



^771^7172^727 rj^^^,r^^^, 



71)72 



(-£^7 ~ -^7') (-^71 ^ -^7') (-^72 ~ ^7') 



(138) 
(139) 
(140) 



(141) 



(142) 

(143) 
(144) 

(145) 
(146) 
(147) 



(148) 



In the end of this subsection, we would like to point out that one of purposes introducing the decompositions of 
S'-product and the contractions and anti-contractions of g-product is to eliminate the apparent (not real) singularities 
and find out all the limitations from the contributions of g'-product contractions. It is important to express the results 
with the physical significance. 
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VII. IMPROVED FORMS OF PERTURBED SOLUTION OF DYNAMICS 



In fact, the final aim using the decomposition and then obtaining contraction and anti-contraction is to include 
partially contributions from the high order approximations of nondiagonal part of Hi in our improved scheme of 
perturbation theory. In this section, making use of the skill to calculate contractions and anti-contractions of g- 
product, we can obtain the zeroth, first, second and third order improved forms of perturbed solutions including 
partial higher order approximations. 

In mathematics, the process to obtain the improved forms of perturbed solutions is a kind of skill to deal with 
an infinity series, that is, according to some principles and the general term form to rearrange those terms with the 
same features together forming a group, then sum all of the terms in such a particular group that they become a 
compact function at a given precision, finally this infinity series is transformed into a new form that directly relates 
with the studied problem. More concretely speaking, since we concern the system developing with time t, we take 
those terms with {—iyit)e~^^'* , (— iyit)^e~'^'*/2' and {—iyit)^e~^^'*/3l, ■■■ with the same factor / together forming 
a group, then sum them to obtain an exponential function /exp[— 1(2:; + yi)i\. The physical reason to do this is 
such an exponential function represents the system evolution in theory and has the obvious physical significance in 
the calculation of transition probability and perturbed energy. Through rearranging and summing, those terms with 
factors f^e"'^"***, (a = 1, 2, • • • ) in the higher order approximation are added to the improved lower approximations, 
we thus can advance the precision, particular, when the evolution time t is longer. We can call it "dynamical 
rearrangement and summation" technology. 



A. Improved form of the zeroth order perturbed solution of dynamics 



Let us start with the zeroth order perturbed solution of dynamics. In the usual perturbation theory, it is well-known 



(149) 



77 



where a^i — ($''' |^(0)). Now, we would like to improve it so that it can include the partial contributions from higher 
order approximations. Actually, we can find that ^2(2) and A^^nn^c) have the terms proportional to (—it) 



(-it)c-'^-* 
(-ii)c-'^^* 



E- 

. 71 

E 



.71 .72 



1 



1^771 |2 



^77' I 



Introduce the notation 



- E 



\9T 



„77l „7l72 „727 

9i ill yi 



71 



E^ E-y 



1 



[E^ - E^J{Ej - E^^) 



„77i „7172„727 

yi di yi ■ 



(150) 
(151) 

(152) 
(153) 



It is clear that G^"-* has energy dimension, and we will see that it can be called the oth revision energy. Let us 
add the terms ((TCT|l and the related terms in A4(te"'^^*, D), A4(i2e"'^^*, D), A5(te-'^^*, D), A^it'^e''^-'^ 

AQ{t^e~^^~i* and >l6(i'^e~'^^*) given in Appendix B together, that is. 



,D), 



477 



til 

2! 



-f(-it) (g(2) + g(3) + gW-i-g(5)) 



(-it) 



2^.(2) (^(4) 

2! ^^7 '-'7 



(154) 



Although we have not finished the more calculations, from the mathematical symmetry and physical concept, we can 
think 



477 



"i-f (-it) (g(2) + g(3) + g(^) + g(5)) 

+^ (Gf + Gf + G(4) + G(5)) ' + . • 



(155) 
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New terms can appear at A^, As, Ag and Aiq. So we obtain the improved form of zeroth order perturbed solution of 

dynamics 



77' 



(156) 



It is clear that clf^ is real. In fact, G'^^' is also real. In order to prove it, we exchange the dummy indexes 71 and 

(3) 

72 and take the complex conjugate of G7 , that is 



71.72 



1 



(^7 ~ -^71) (-^7 ~ -^72) 



„77i „7172 _727 

yi 9i Ui 



(157) 



where we have used the relations = gf^^^ for any (3i and /?2 since Hi is Hermit. Similar analyses can be 

applied to and . These mean that e \ y~^y~^y~^'r)is still an oscillatory factor. 

B. Improved form of the first order perturbed solution of dynamics 

Furthermore, in order to include the partial contributions from higher than zeroth order approximation, we need 
to consider the contributions from nondiagonal elements in the higher order approximations. 
Well-known usual first order perturbed part of dynamics is 



= E 



Q—iE^t 



^7"'l^^) = E 



7,7' 



77 

Si 



|$T). (158) 



It must be emphasized that H\ is taken as only with the nondiagonal part for simplicity. That is, we have used the 
skill one stated above. 

Thus, from A^{t%-'^'^\^) and Ai{tQ-'^'^\^), Ai{ee-'^'^\Vi), A^{t^e-'^'^\W), A(,{t'Q-'^-^\^) in the Appendix 
B, it follows that 



All 
^1:1 



{E^ - Ey) 



1 



+ 



-i-it) (G(2) + Gf +GW) 
%^2GfGf + ••• 

+ i^it) {g^} + G^ + d^) 



+ 



{Ej — Ey 

(-ii)! M2)\ 2 (-ii)" 0^(2)^(3) 



2! 



+ 



2, 2G^1^G^1^ + ... 



(159) 



Therefore, the improved perturbed solution of dynamics is just 

= E 

7j7' 



7' 



— E' 



1' 



77 

gV 



(160) 
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C. Improved second order- and third order perturbed solution 

Likewise, it is not difficult to obtain 



■E 

7)717' 



+ 



+ 



„77i„7i7': 

9i 9i 0. 



77 



(£//y — -^71 ) (-^7 — -£/7' ) (-^7 — -^71 ) (-^71 — -E'7' ) 
57''^ffr»?77'>«7'l*^)- 



(-2^7 ~ -^7') (-^71 ~ -^7') 



i(B^+G^^>)t 



-i(B^+G^^')t 



+ 



[E^ - E^i ) {E^ - E^^ ) (E^ - E^i ) - Ej^ ) 



(-5^7 - -^71 )^ (-^71 ~ -^72 ) (-^7 ~ ^l^f {Eji - E^2 ) 



7>7 V 71 L 



eTe 

7.7' 

E 



i(s^+G^^))t 



-i(B,+G(,=')t 



{E^ — -^71) (-^7 — -^7') (-^7 — -^71) (-^7 — -^7') 



„77l „7l7„77 
i/l i/1 i/l 



^772 



-i(£;^,+G^=^^))t 



??7i 



(-S7 - -E^i) (-B^ - E^^) {Ej - Ey) {E^ - E^^) {Ej^ - Ej^) {E^^ - E^>) 



+ 



^772 



{E^ - E^^) [Ery^ - Ej^) {Ej^ - Ey) 



(161) 



(162) 



D. Summciry 

Obviously, our improved form of perturbed solution of dynamics in the third order approximation is 

|*(t)) = ^|vI/W(t))j + 0(ij4). (163) 

i=0 

However, this solution including the contributions from the whole AJ'' {te), Al^ (i^e) parts up to the fifth order 

approximation and the whole JCp (t^e), JCp (t^e) parts in the sixth order approximation. After considering the 
contractions and anti-contractions, we see the result corresponds to the replacement 

in the JCp (e) part, where 

^7. =^7, +E^7"^' (165) 

a=2 

i = 0,1,2,---, and 70 = 7- Although the upper bound of summation index a is different from the approximation 
order in the finished calculations, we can conjecture it may be taken to at least 5 based on the consideration from the 
physical concept and mathematical symmetry. For a > 5, their forms should be similar. From our point of view, such 
form is so delicate that its form happens impossibly by accident. Perhaps, there is a fundamental formula within it. 
Nevertheless, we have no idea of how to prove it strictly and generally at present. 

Actually, as soon as we carry out further calculations, we can include the contributions from higher order approx- 
imations. Moreover, these calculations are not difficult and are programmable since we only need to calculate the 
limitation and summation. Therefore, the advantages of our solution have been made clear in our improved forms of 
perturbed solution of dynamics above. In other words, they offer clear evidences to show it is better than the existing 
method in the precision and efficiency. In the following several sections, we will clearly demonstrate these problems. 



(164) 
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VIII. IMPROVED TRANSITION PROBABILITY AND REVISED FERMI'S GOLDEN RULE 



One of the interesting applications of our solution is the calculation of transition probability in a general time- 

indcpcndcnt quantum system. It improves the well-known conclusion because our solution include the contributions 
from the high order approximation. Moreover, in terms of our improved forms of perturbed solution, it is easy to 
obtain the high order transition probability. In addition, for the case of sudden perturbation, our scheme is also 
suitable. Since considering time-independent system, it is not necessary to using the interaction picture. 
Let us start with the following perturbed expansion of state evolution with time t, 



(166) 



1=0 7 



When we take the initial state as |^''^), from our improved first order perturbed solution, we immediately obtain 



(167) 



where 



(168) 



Here, we use the subscript "I" for distinguishing it and the usual result. Omitting a unimportant phase factor e ''^t*, 
we can rewrite it as 



T'l K - E, 



^1 _ e'^^''*^ 



(169) 



where 



uj^l3 = E.y — Ep. 



(170) 



Obviously it is different from the well known conclusion 



.(1) 



(1 _ e''^7/3t) ^ 



(171) 



where 



^ifs = — Ep. 



(172) 



Therefore, our result contains the partial contributions from the high order approximation. 
Considering the transition probability from |$^) to 1$''') after time T, we have 



yd 



7/3 



1 _ e''^-v/3T 



7/3 

9i 



2 sin^ (t^7/3T/2) 
K/3/2)' 



(173) 



In terms of the relation 



sin^ X — sin^ ^ ~ 2 1^°^^^^) ~ cos(2a;)] , 



(174) 



we have the revision part of transition probability 



2 cos {cOj/sT) — COS {cOjfjT) 



(175) 



If plotting 



s in^ _ fuj^py sin^ 



K/3/2)' 



V'^7/3/ (w^/3/2) 



(176) 
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we can see that it has a wen-defined peak centered at oj^p = 0. Just as what has been done in the usual case, we can 
extend the integral range as — oo — *■ oo. Thus, the revised Fermi's golden rule 



W = ?i'F+ Aw, 



where the usual Fermi's golden rule is 



10 
9i 



(177) 



(178) 



in which w means the transition velocity, p{E-y) is the density of final state and we have used the integral formula 



■ 2 

sm X 



while the revision part is 



Aw = 2 dE^p{E^) 



7/3 

9{ 



cos (uj^pT) — cos {uj^fiT) 



(179) 



(180) 



It is clear that uj-yp is a function of S-y, and then a function of w-y^. For simplicity, we only consider uj^p to its second 
order approximation, that is 



71 



\9T 



9i 



0{Hl). 



(181) 



Again based on dE^ ~ dw-^^, we have 

/•oo 

Aw = 2 duj-ypp {lo^p + Ep) 



7/3 

9i 



2 cos [^7/3^] — cos [w {ijJ^p) T] 

T{uj^pf 



(182) 



It seems it is not easy to deduce the general form of this integral. In order to simplify it, we can use the fact that 
uj^p — LUryp is a smaller quantity since 



AUJjP = LO-yp — UJ-yp 



1=2 



For example, we can approximatively take 

cos {ujfpT) - cos (oj-ypT) « T {oj^p - uj^p) sin (uj^pT - uj^pT) , 



(183) 



(184) 



then calculate the integral. We will study it in our other manuscript (in preparing). 

Obviously, the revision comes from the contributions of high order approximations. The physical effect resulted 
from our solution, whether is important or unimportant, should be investigated in some concrete quantum systems. 
We will discuss them in our future manuscripts (in preparing). 

It is clear that the relevant results can be obtained from the usual results via replacing uj^p in the exponential 
power by using uj-yp. To save the space, we do not intend to mention them here. 

In fact, there is not any difficult to obtain the second- and three order transition probability in terms of our improved 
forms of perturbed solution in the previous section. More higher order transition probability can be given effectively 
and accurately by our scheme. 

IX. IMPROVED FORMS OF PERTURBED ENERGY AND PERTURBED STATE 

Now we study how to calculate the improved forms of perturbed energy and perturbed state. For simplicity, we 
only study them concerning the improved second order approximation. Based on the experience from the skill one 
in Sec. IVII we can, in fact, set a new E and then use the technology in the usual perturbative theory. That is, we 
denote 



E^ 



E^, + Gi^.^ + G^^l 



(185) 
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7,7' ^. 



E 

71 



71 



771 717 r 
2 i/l Hi "77 



(-E^ — ) (-£'7 — -E^7') (-^7 ^ -^71) (-^71 ~ -^7') 



57"^5r'777' «7'l'i>") + 0(i?f)- 



Because that 



we have 



|vI'(t))j = ^e-'^-*<57ya7H$^), 
7.7' 



(186) 



(187) 



i?7 — i?7 



E^ -~ E-« 



71 



-57^' "E- 

71 {E^-EjiY 



^7 -^71 „77l „7l7r 

TP •^'i / -TP N2i'l yi "77' 

ii7 — Hj^- '-^ ' '■ * 



{E-y — -B71 ) (-^7 — -^7') (-^7 ^ -^71) (-^71 ~ -^7') 



771 717 I 
W 5l '777' r °^7'- 



(i?7 — -E^7') (-B71 — ^7') 

In the usual perturbation theory, Hi is taken as a perturbed part with the form 

Hi ~ \Hi, 



(188) 



(189) 



where A is a real number that is called the perturbation parameter. It must be emphasized that E^. can be taken as 
explicitly independent perturbed parameter A, because we introduce A as a formal multiplier after redefinition. This 
way has been seen in our skill one. Without loss of generality, we further take Hi only with the nondiagonal form, 
that is 



^7172 ^ ^^7172 ^XgJ'^\ 



(190) 



Then, we expand both the desired expansion coefficients and the energy eigenvalues Et in a power series of 
perturbation parameter A: 



Er ^ E^'4'l 



1=0 



(191) 
(192) 



1=0 



1. Improved 0th approximation 

If we set A = 0, ea. l(TFS|) yields 

p(0)„(0) _ p „(0) (IQ'}\ 

where 7 runs over all levels. Actually, let us focus on the level 7 = /?, then 

4°l = Ep- (194) 

When the initial state is taken as \^^), 

4a = ^7/3- (195) 

Obviously, the improved form of perturbed energy is different from the results in the usual perturbative theory because 
it includes the contributions from the higher order approximation. However, the so-call improved form of perturbed 
state is the same as the usual result. 
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2. Improved 1st approximation 

Again from eq. l(TS5jl it follows that 



p(0)„(l) , p(l)„(0) _ p „(1) , £'7 - -5'7' ~77' (0) nofiA 



^7 -^7' 
7 ' ' 



When 7 = /3, it is easy to obtain 



If 7 ^ /3, then 



= 0. (197) 



aW = _^ 1 gi/. (198) 

It is clear that the first order results are the same as the one in the usual perturbative theory. 



S. Improved 2nd approximation 



Likewise, the following equation 



p(2)„(0) „(1) , p(0)„(2) _ p „(2) , -E7 --g7' ~77' (1) 



E 



7' 



71 ^ 



71,7 



E^ 



{E^ — £71) (-£71 ^ -£7') (£7 ^ -£7') (£71 — £7') 



771 717 (0) 
5i 5i %7'«7';i- 



(199) 



is obtained and it yields 



E, 



(2) 
'T;I 



0, 



(200) 



if we take 7 = /3. When 7 ^ /3, we have 



(2) 

^7;i = 



E 

71 



(E^ - Ep) (E^i - E^) 



^7l~7l/3 

5l ffl ??7/3- 



(201) 



It is consistent with the nondiagonal part of usual result. In fact, since we have taken H^'^ to be nondiagonal, it 
does not have a diagonal part. However, we think its form is more appropriate. In addition, we do not consider the 
revision part introduced by normalization. While E^^.j = is a new result. 



4- Summary 

Now we can see, up to the improved second order approximation: 

Et,/5 ^Efj^Ep + Gf + Gf . (202) 

Comparing with the usual, they are consistent at the former two orders. It is not strange since the physical law is 
the same. However, our improved form of perturbed energy contains a third order term. In other words, our solution 
might be effective in order to obtain the contribution from high order approximation. The possible physical reason is 
that a redefined form of the solution is obtained. 
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In special, when we allow the H'l'^ to have the diagonal elements, the improved second order approximation becomes 

Ets ^Ep + h1 + Gf + Gf . (203) 

Likewise, if we redefine 

E^, = E^, + G^2) + Gi3) + . (204) 
Thus, only consider the first order approximation, we can obtain 

Et,0 ^Ep + hf + Gf + Gf + Gf . (205) 
In fact, the reason is our conjecture in the previous section. The correct form of redefined E^. should be 

Et..p « + hf + Gf + Gf + Gf +Gf+---. (206) 

This implies that our improved scheme includes the partial even whole significant contributions from the high order 
approximations. In addition, based on the fact that the improved second approximation is actually zero, it is possible 
that this implies our solution will fade down more rapidly than the solution in the usual perturbative theory. 

Actually, the main advantage of our solution is in dynamical development. The contributions from the high order 
approximation play more important roles in the relevant physical problems such as the entanglement dynamics and 
decoherence process. For the improved perturbed energy, its high order part has obvious physical meaning. But, for 
the improved form of perturbed state, we find them to be the same. 



X. EXAMPLE AND APPLICATION 



In order to concretely illustrate that our exact solution and perturbative scheme are indeed more effective and more 
accurate, let us study an elementary example: two state system, which appears in the most of quantum mechanics 
textbooks. Its Hamiltonian can be written as 

where we have used the the basis formed by the unperturbed energy eigenvectors, that is 

l*^)=(;). I*^)=(?)- (208) 

In other words: 

Ho\<l>-') = E^\<l>->), (7 = 1,2) (209) 

where 

""(I'D- w 



Thus, this means the perturbed Hamiltonian is taken as 



This two state system has the following eigen equation 

H\^^) = E^\^^). (212) 
It is easy to obtain its solution: corresponding eigenvectors and eigenvalues 

1*2) = ^ ( ^^21 - '^2^1 ^ ; (214) 
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= l{E, + E2+co^^); 



(215) 

(216) 



where = |Vi2| = |V2i|, LO21 = E2 — Ei, lo^i = -Ej — Ef = ^/4\V^'+7o^, and we have set E2 > Ei without loss of 
generality. 

Obviously the transition probability from state 1 to state 2 is 

2 



P'^{1 ^2) = |($ 

= \v\ 



$1)1 



71,72 = 1 



, sin^ Kit/2) 
Ki/2)2 ■ 



(217) 



In the usual perturbation theory, up to the second order approximations, the well-known the perturbed energies 
are 

\V\^ 



El 



E^ = Ei + 



W21 
W21 ' 



While, under the first order approximation, the transition probabihty from state 1 to state 2 is 
Using our scheme, only to the first approximation, the corresponding the perturbed energies are 



El = El 
E2 = El + 



W21 

1^1! ^ 



^21 

\v\^ 



where we have used the facts that 



(2) _ JVl^, 
0J21 
- = 0, 

4 



-g: 



(2) 



2 ) 



2 

(4) _ \V\- _ _ (4) 



3 - ~<^2 

'21 



(218) 
(219) 

(220) 

(221) 
(222) 

(223) 
(224) 
(225) 



Obviously, under the first order approximation, our scheme yields the transition probability from state 1 to state 2 
as 

2 , 



Pi(l -^2) = \V 



2 sin {uj2it/2) 



(u;2i/2)2 • 

where a;2i = E2 — Ei. Therefore we can say our scheme is more effective. Moreover, we notice that 



El = El — 



+ 



U!21 OJ21 

= Ei + 0{\V\^) 
\V\' 



+o{\vn 



= Er+ 



+oi\vn 



■'21 



E2 = £;i + ^-^ + o(|yH 
W21 W21 

= E2 + 0{\Vf) 

11/14 



E. 



2 , ,3 
'^21 



+o{\vn 



(226) 

(227) 
(228) 
(229) 

(230) 
(231) 
(232) 
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and 



P'^(l ^ 2) = \V\ 



I sin^ (w2ii/2) 

(t^2l/2)2 



sin (cL'2it) sin^ (a;2iV2) 
2(w2i/2)2 " (u;2i/2)3 



(W21 - W2l) + 0[(W21 - W2l) ] 



Pi{l^2)-\V\- 



, sin^ (tj2ii/2) 
(a;2i/2)3 



(tJ21 - ^^2l) + C[(^^21 - ^21)^ 



P(l 



sin(a;2i0 sir? (u}2it/2) 



_2(c^2i/2)2 (c^2i/2)3 
Therefore, we can say that our scheme is more accurate. 



(aj2i - W21) + 0[{uj2i — ^21) 



(233) 
(234) 
(235) 



XI. CONCLUSION AND DISCUSSION 



In the end, we would Uke to point out that our general and explicit solution l|35|l or its particular forms (|42I43|I of 
the Schrodinger equation in a general time-independent quantum system is clear and exact in form in spite of it is an 
infinity series. It must be emphasized that its distinguished feature is to first express the exact solution in a general 
time-independent quantum system by using the c-number matrix elements rather than an operator form. We can say 
that our exact solution is more explicit than the usual nonperturbed solution of the Schrodinger equation. Moreover, 
our deducing methods give up some preconditions used in the usual scheme. So we can say our solution is more 
general. Just as its explicit and general form, our exact solution not only has the mathematical delicateness, but also 
can contain more physical content, obtain the more efficiency and higher precision and result in new applications. 

We obtain our exact solution based on our expansion formula of power of operator binomials and the matrix 
representation of time evolution operator, which is a new way to study the quantum dynamics and the perturbation 
theory. Its idea is different from the usual one. At our knowledge, this formula is first proposed and strictly proved. 
Besides its theoretical value in mathematics, we are sure it is interesting and important for expressing some useful 
operator formula in quantum physics. 

In the process of deducing our exact solution, we prove an identity of fraction function. It should be interesting in 
mathematics. Perhaps, it has other applications to be expected finding. 

By virtu of the idea of quantum field theory, we build the concepts of contraction and anti-contraction using in the 
limitation computation. Moreover we develop two skills to deal with an infinity series, in special, so-called dynamical 
rearrangement and summation. These methods have played important roles in our scheme. We believe they can be 
useful technologies in mathematical and physical calculation. 

As an example, we give out the concrete form of solution when the solvable part of Hamiltonian is taken as the 
kinetic energy term in a general quantum system. Its aim is to account for the fact that our solution perhaps extends 
the applicable range of perturbation theory and lose its preconditions. 

We can see that in our solution, every expanding coefficient (amplitude) before the basis vector 1$''') has some 
closed time evolution factors with exponential form, and includes all order approximations so that we can clearly 
understand the dynamical behavior of quantum systems. Although our solution is obtained in representation, its 
form in the other representation can be given by the representation transformation and the development factors with 
time t in the expanding coefficients (amplitude) do not change. 

Because we obtain the general term form of time evolution of quantum state, it provides the probability considering 
the partial contributions from the high order even all of order approximations. After developing the contraction and 
anti-contraction skills and dynamical rearrangement and summation technology, we realize this probability by present 
a perturbative scheme. From our exact solution transferring to our perturbative scheme is physically reasonable and 
mathematically clear. This provides the guarantee achieving high efficiency and high precision. Through finding the 
improved forms of perturbed solutions of dynamics, we generally demonstrate this conclusion. Furthermore, we prove 
the correctness of this conclusion via calculating the improved form of transition probability, perturbed energy and 
perturbed state. Specially, we obtain the revised Fermi's golden rule. Moreover, we illustrate the advantages of our 
exact solution and perturbative scheme in an easy understanding example of two state system. All of this implies 
the physical reasons and evidences why our exact solution and perturbative scheme are formally explicit, actually 
calculable, operationally efficient, conclusively more accurate. 

From the features of our solution, we believe that it will have interesting applications in the calculation of en- 
tanglement dynamics and decoherence process as well as the other physical quantities dependent on the expanding 
coefficients. Of course, our solution is an exact one, its advantages and features can not be fully revealed only via the 
perturbative method. 

Undoubtedly, the formalization of physical theory often has its highly mathematical focus, but this can not cover 
its real motivations, potential applications and related conclusions in physics. Our exact solution and perturbative 
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scheme are just so. In order to account for what is more in our solution and reveal the relations and differences 
between our solution and the existing method, we compare our solution with the usual perturbation theory. We find 
their consistency and relations. In fact, our solution has finished the task to calculate the expanding coefficients of 
final state in Hq representation and obtain the general term up to any order using our own method. But the usual 
perturbation theory only carries out this task from some given order approximation to the next order approximation 
step by step. In a sentence, more explicit and general feature of our exact solution can lead to more physical conclusions 
and applications. 

It is worth pointing out that our solution, different from the time-dependent perturbation theory, presents the 
explicit solution of recurrence equation of expanding coefficients of final state in Hq representation, but we pay for 
the price that H is not explicitly dependent on the time. In short, there is gain and there is lose. 

In fact, a given lower order approximation of improved perturbation solution including the partial contributions 
from the higher order even all of order approximations is obtained by rearranging and summing a kind of terms, just 
like it has been done in the quantum field theory for the particular contributions over a given type of the Feynman's 
diagrams. It is emphasized that these contributions have to be significant in physics. Considering time development 
form is our physical idea and adding the high order approximations with the factors t"e~'^^i*, (a = 1, 2, • • • ) to the 
improved lower order approximations definitely can advance the precision. Therefore, our dynamical rearrangement 
and summation technology is appropriate and reasonable from our view. 

For a concrete example, except for some technological and calculational works, it needs the extensive physical 
background knowledge to account for the significance of related results. That is, since the difference of the related 
conclusions between our solution and the usual perturbation theory is in high order approximation parts, we have to 
study the revisions (difference) to find out whether they are important or unimportant to the studied problem. In 
addition, our conjecture about the perturbed energy is based on physical concept and mathematical consideration, it 
is still open at the strict sense. As to the degenerate case, except for the complicated expressions, there is no more 
new idea. 

Based on the above statements, our results can be thought of as theoretical developments of quantum dynamics. 
The extension of our solution to mixed states is direct, and the extension of our solution to open systems see our 
manuscript 

It must be emphasized that the study on the time evolution operator plays a central role in quantum dynamics. 
From our point of view, one of the most main results in our method is to obtain a general term form of any order 
of Hi (perturbative part of Hamiltonian) for the time evolution operation in the representation of Hq (solvable part 
of Hamiltonian), however, the usual perturbation theory has not really finished it and only has an expression in the 
operator form. Because the universal significance of our new expression of time evolution operator, we wish that it 
will have more applications in quantum theory. Besides the above studies through the perturbative method, it is more 
interesting to apply our exact solution to the formalization study of quantum dynamics in order to further and more 
powerfully show the advantages of our exact solution. 

In summary, we can say our present results are helpful to understand the theory of quantum dynamics and provide 
some powerful tools in the calculation of quantum dynamics. More investigations are on progressing. 
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In this appendix, we would like to prove our identity (|26|l . For simplicity in notation and universality, we replace 
the variables E^^ by Xi as well as E[^,l] by x[l]. It is clear that the common denominator D{x[n]) in the above 
expression (|26|l (the index I is replaced by n) reads 
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APPENDIX A: THE PROOF OF OUR IDENTITY 



n n 



D{x[n])^l[ He 



(Al) 



i=l j=i 



while the i-th numerator is 



ni{x[n],K) 



di{x[n\) 



.K 

■i ; 



(A2) 
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and the total numerator N{x[n], K) is 



n+l 



ni{x[n],K). 



(A3) 



In order to simplify our notation, we denote ni(x[n]) — ni(x[n], 0). Again, introducing a new vector 

= {x2,X3, • • ■ ,a;„,a;„+i}, 
= {xi,--- ,a;i_i,a;i+i,--- ,a;„+i}, 
^n+iW] = {a;i,X2, • • • ,a:„_i,a;„}. 



(A4) 
(A5) 
(A6) 



Obviously, xf{x[n]) with n components is obtained by deleting the i-th. component from x\n]. From the definition of 

ni{x[rL\), it follows that 



n,{x[n])^D{xY[n]) 



(A7) 



Without loss of generality, for an arbitrary given i, we always can rewrite = y[n — 1] = {yi, j/2, • • • , J/n-i, Un] 

and consider the general expression of D{y[n ~ 1]). It is easy to verify that 



Diy[l])^-y2n,{y[M)+yiMy[M), 
D{y[2])^y2y3ni{y[2]) - yiy3n2{y[2]) + yiy2n3{y[2]). 

Thus, by mathematical induction, we first assume that for n > 2, 

n 

D{y[n - 1]) = Y.{-lY'~'^+^"-'^pMn ~ - 1]), 

4=1 

where we have defined 



Vj- 



(A8) 
(A9) 



(AlO) 



(All) 



As above, we have verified that the expression l|ITO)l is valid for n = 2, 3. Then, we need to prove that the following 
expression 



n+l 



D{v[n]) = ^(-l)(*-i)+>.(2/M)n,(2/[n]) 



(A12) 



i=l 



is correct. 

To our purpose, we start from the proof of a conclusion of the precondition HA10() as the following: 



n+l 



(A13) 



According to the relation (|A7|) and substituting the precondition (|A10|) . we have 



n+l 



n+l 



Y^{-iy-'nMn]) - 5](-l)*-ii^(2yf [r 



i=l 
n+l 



i=l 



J2i-l)^'-'^^^''-'^MyF[n])n,{yF[n]) 



n n+l 



(A14) 
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Based on the definitions of pi{x[n]) and ni{x[n]), we find that 



MyF[n])^\ ^ ^ ^ , (A15) 

rn,_i(yf[n]) (If z > j) 

^.(yfN)= , , , ■ (A16) 

in,(yf+Jn]) (Ifz<j) 

Thus, the right side of eQ. (|A14|) becomes 

n n+l 

j=i i=i 

n n+1 

= (-i)""'E E (-i)"+^p.(2/fN)«.(yfN) 

n j 
n n+1 

= (-i)"-'E E (-i)'+^P.-i(yfN)«.-i(yfW) 

+(-i)"^'EE(-i)'^'p^(^f+iN)«.(yf+iN)- (A17) 

Setting that i — 1 ^ i for the first term and j + 1 ^ j for the second term in the right side of the above equation, we 
obtain 

(-1)""' E E(-i)^^'^'^(2^«''N)«.(yf N) 

n n 

= (-i)"-'EE(-i)'^'"V.(yfW)n.(2/fW) 
EE(-i)'^'"'^'«(2^fN)".(yfN) 

j=2 ^=l 

n n 

= (-i)"EE(-i)"^'^''(yfW)"«(2^fW) 

n+l 

+ (-!)" ^(-l)''+("+i)K(y^+i[n])n,(y^+i[n]). (A18) 
Again, setting that i ^ j in the right side of the above equation gives 

n n+l 

(-i)""'EE(-i)*^'p^(2^«''w)"^(2^»''N) 

= (-i)"EE(-i)'^'p^-(y^''W)'^^(2^«''N)- (A19) 

Thus, it imphes that 

n n+l 

(-l)""'EE(-l)'^'^'^(^f N)n,(yf [n]) = 0. (A20) 
j=i i=i 
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From the relation (|A14|I it follows that the identity ljA13|) is valid. 

Now, let us back to the proof of the expression (|A12() . Since the definition of D{y[n]), we can rewrite it as 

D{y[n])^D{y[n~l])d„+Mn]). 
Substituting our precondition (jA10|l yields 

n 

D{y[n]) - - - 

i=l 

Note that 

nt{y[n - l])dn+i{y[n]) = n,{y[n]){y, - y^+i), 

Pt{y[n - l])yn+i = Pt{y[n]), 
Pi{y[n-i])yi = Pn+i{y[n]), 

we have that 

n n 
n r n ~ 

= Y.'^-lY'''^+y^{y[n])nMn]) + (-1)- v+i(y w) E(-i)^"'"^(y W) 

n 

= E(-l)*'"'^^>^(yW)^»(j/W) + (-1)"- V+i(2/W) [-(-l)"n„+i(y[n])] . 
1=1 

In the last equality we have used the conclusion (jAlSp of our precondition, that is 

n 

J2{-iy-Wy[n]) = -(-l)"n„+i(yN). 

Thus, ea. HA26|l becomes 

D{y[n]) = J2{-lf-'^+"pMn])n^{y[n]). 



4=1 



(A21) 



(A22) 



(A23) 
(A24) 
(A25) 



(A26) 



(A27) 



n+l 



(A28) 



The needed expression (jA12ll is proved by mathematical induction. That is, we have proved that for any n > 1, the 
expression (|A12|) is valid. 

Since our proof of the conclusion (|A13I) of our precondition is independent of n (n > 1), we can, in the same way, 
prove that the identity (|A13p is correct for any n > 1. 

Now, let us prove our identity (|26|l . Obviously, when K =: we have 



ji+i 



Yi-iy-' — - 



1 



d,{x[n]) D{x[n]) 



Vi+l 



= 0, 



(A29) 



where we have used the fact the identity (|A13|I is valid for any n > 1. Furthermore, we extend the definition domain 
of {x[n]) from K > n to K > Q, and still write its form as 



n+l 



di{x[n\) 



Obviously, ea. (jA29|) means 



C°(x[n]) = 0, (n > 1). 



(A30) 



(A31) 
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In order to consider the cases when i^T 7^ 0, by using of di{x[n]){xi — Xn+2) = di{x[n + 1]) {i < n + 1), we obtain 



di{x[n]) 



1=1 

"+1 ^a: 



d,(a;[n+ 1]) 
-x„+2^(-l)'-i 



n+l 7^ 



.^^ d,(a;[n + l]) 

It follows the recurrence equation as the following 

C^ix[n]) = C^+,'i4n + 1]) - x„+2C,f+i(x[n + 1]). (ASS) 

It implies that since C^{x[n]) = for any n > 1, then Cj^^i{x[n + 1]) = for any n > 1 or C^{x[n]) ~ for any 
n > 2; since C}i{x[n]) = for any n > 2, then C^^_f_i{x[n + 1]) = for any n > 2 or C^{x[n]) = for any n > 3; ■ ■ ■ , 
since C^{x[n]) = for any n > {k + I) {k > 0), then C';+l{x[n + 1]) for any n > (fc + 1) or C!^+\x[n]) = for 
any n > (fc + 2); • • ■ . In fact, the mathematical induction tells us this result. Obviously 

C^{x[n]) = 0, {UO<K<n). (A34) 

Taking K = n in ea. (|A33|) and using ea. (|A34|) . we have 

C:{x[n])^C:+l{x[n + l])^l, (n > 1), (ASS) 

where we have used the fact that Cj'^(x[l]) — 1. Therefore, the proof of our identity H2t)|) is finally finished. 



APPENDIX B: THE CALCULATIONS OF THE HIGH ORDER TERMS 



Since we have taken the Hi only with the non diagonal part, it is enough to calculate the contributions from them. 
In Sec. IVII the contributions from the first, second and third order approximations. In this appendix, we would like 
to find the contributions from the fourth to the sixth order approximations. The used calculation technologies are 
mainly to find the limitation, dummy index changing and summation, as well as the replacement gl''^' Tj^i-y^ = gl''^^ 
since gj^'^' has been nondiagonal. They are not difficult, but are a little lengthy. 



1. 1 = 4 case 



For the fourth order approximation, its contributions from the first decompositions consists of eight terms: 

A^' = A2^' (ccc) + AJ^' (ccn) + AJ^' (cnc) + AJ*' {ncc) 

+A'P' [cnn) + AJ^' {ncn) + AJ^' [nnc) + AJ^' [nnn). (Bl) 

We will see that the former four terms have no the nontrivial second contractions, the fifth and seven terms have one 
nontrivial second contraction, 

A'P' {cnn) = AJ^' {cnn,kc) + A'P' {cnn\kn), (B2) 
A'P' {ncn) = AJ^' {ncn,c) + AJ^' {ncn,n), (BS) 
A'P' {nnc) = AJ^' {nnc, ck) + AJ^' {nnc, nk). (B4) 
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In addition, the last term in ea. (|Bl|) has two nontrivial second contractions, thus 

/ / / / / 

Ap (nnn) — A'P {nnn,cc) + Ap {nnn,cn)A2'^ {nnn,nc) + Ap {nnn,nn), 

its fourth term has also the third contraction 

AP {nnn, nn) = AP {nnn, nn, c) + Ap {nnn, nn, n). 



(B5) 



(B6) 



All together, we have the fifteen expressions of the contributions from whole contractions in the fourth order approx- 
imation. 

First, let us calculate the former four terms only with the first contractions and anti-contractions, that is, with 
more than two 6 functions (or less than two rj functions) 



AP (ccc) 



E 

71, •■ 

E 



-iE^.t 



71 



d^{Eh,4])_ 
3e-*^^i* 



9i 



'^7fc7fc + 2 ) '^7l7'^757' 



\fc=l 
2Q-iE-,t 



{E^ - E^J^ {Ej - EjJ* ^ \ej - E^-^f 



-(-it) 



{-uy 



(^7~^7i) ^ (£^7^-^71)^ 



(B7) 



Ap'{ccn) = E(-l)'~' 

7l.--- ,75 Li=l 

= E 



Q — iE-y^t 



4 

n 



„7j7j + l 

9i 



^71 73 ^12 74 ^773 75 '^71 7 "^75 7' 



(i?-y — -^^7i) (-^7 — -^7') (-^7 — -^71 ) (-^7 — ^7') 



{E^ - Eyj'^ {E^^ - Ey)'^ {Ey ~ i?7i)^ (£'71 - i?7') 



^7 -^71/ V-^71 -^7 

(-E-y — -E7') (-£^71 ~ -E^7') 



'7 ^71 / V-^71 -^7 



(-E^ — E^^) (£71 ~ ^7' 



(-E^ — ) {Ery — E^l 



(B8) 



AP (c 



E 

7l,-- - ,75 



E 



E(-i)'-^ 



rf.(i?[7,4]) 



■J~J „7j7j + l 



9i 



i=i 



'^71 73 74 '^73 75 ^71 7 "^75 7' 



(E^ - E-yi ) (E^ - E-y2 )^ (E-y - E^if {E-i - E^^f 



71 ,72 L ^ Ti ^ T ''s ^ 

(E^ - E^ {E^ - E^2 ) (E-y - E^^ )^ (E^j - E^^ ) 

{-it) 



{Ef - EyJ (E^^ - E^J 

-(-it) 



(E^ - Ey-^ ) (E^ - E^2 ) 
(-it)2 e-'^^* 



(E^ - E^J' (E^ - E^J ' 2 (E^ - E^ J (E^ - E^J 

X 1.97''' 1^ 157''' I^'77l72'577'- 



(B9) 
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AP {ncc) 



E 

71, •• 

E 



E(-i)'-^ 



-lE^-t 



n 



„7j7j + l 

yi 



71 



+ 



^71 73 '^72 74 "^73 75 '^71 7 '^75 7' 



— E'^^^ (-^7 — -^7')^ (-^7 — -^71 ) (-^71 — -^7')^ 



71 V-^7 -^7 



e '-^Ti 



2e- 



(^/'y — -£'71 ) (-^71 — -^7') (-^7 — -^7') (-^71 — -^7') 

(i/^ — -^7') {-^71 — -^7') {-^7 — -^71) (-^71 — -^7') 

2 



-(-«) 



e T 



(iJ-y — -E-y') (£'71 — -Ey' 



^717 



57^^ 5?^ ^77' 



(BIO) 



Then, we calculate the three terms with the single first contraction, that is, with one 5 function. Because one 5 
function can not eliminate the whole apparent singularity, we also need to find out the nontrivial second contraction- 
and/or anti-contraction terms. 



7i,--- >75 



[cnn, kc) = 

1," 



.i=l 



2e~^-^7* 



'^71 73 '772 74 '773 75 '^72 7 ' '^71 7 5 7 



{E^ — -^71) (-^7 — -^7')^ (-^7 — ■^71)'^ (-^7 — -^7')"^ 



-I- 



e 7' 



(-£^7 ~ -^71)^ (-^71 ~ -^7')^ (-5^7 ~ -^7')^ (-£^71 ~ -^7')^ 



2c 



{E^ — E^f ) (-^71 — -£/7' ) 

-i-it) 



-iE-,t 



{Ej — -£^71) (-£^7 ~ -^7') 



(-E7 — -£'7') (-£71 — E'j') 



(Bll) 



(cnn, A;n) = 

IT" 

= E 



71. ■••.75 Li=l 



E(-i)^-^ 



J|^7.7..l 



•^71 73 ^72 74 ^73 75 '?72 7' ^71 7 ^7 



(£^7 - E^^) (£"7 - £^72) (i?7 - -E7') 





{E^ 


- E^^) (£7 - -£72)^ (£'7 - 

Q—iEyt 


Ej') 


{E^ 


- £71)^ {Ej - £^2) {Ej - 


Ej') 


{E^ 


- E^^ )^ (£71 - £72 ) (£71 

g— iE-,2* 


-Ey) 


{E^ 


- E^^f (-2^71 ~ -^72) (-^72 
g-iByt 


-Ey) 


{E^ 


- £.y/)^ (£^,1 - Er^,) (£^2 


-Ey) 



+ 



+(-ii) 



-\E.,t 



(E-y - £71 ) (£7 - £72 ) (£7 - Ey) 



^ Ifll'*^'"! gj^^gj^^ '?7i 72^771 7' '?77'- 



(B12) 



40 



AJ'^ {ncn, c) 



E 

7l.--- .• 

E 

7ii72 



-iE^.t 



n 



91 



Q — iEyt 



'?7l73 ^7274 ^773 75 '^71 7 "^757' '^77' 



{Ej - E^^f {Ej - Ej^f {Ej - Ej^ f {E^ - E^^ ] 



(-^7 ~ -^71)^ (-^71 ~ -^72)^ (-^7 ~ -^71)^ (-^71 ^ -^72) 



+ 



(^7 ^ ^72)^ (^71 ~ -^72)^ 



{Ej - Ej^) [E^^ - E^^) 



{Ey - E^^ ) {Ej - E^2 ) 



(B13) 



{ncn,n) = 



7l.--- ,75 



E(-ir^ 



-iB^.t 



E 



Q-iE^t 



4 

n 



„7j73 + i 



^71 73 '^72 74 ^73 75 '^71 7 '^75 7' ^777' 



71,72 ~ -^Ti) (-^7 ~ -^72) (-^7 ~ Ey) 



+ 



+ 



{E, 


- Ej^ ) {E^^ 
e" 


-iE-y^ t 


-Eyf 


{E^ 


- Ej^ ) (iJ^i 
e" 


- -E'72)^ (^'71 




{E^ 


- Ejj'^ {Ej^ 
e" 


- Ej2 ) (i;^! 

"'^72* 


-Ey) 


(^7 


- £^,3) {E^^ 
e" 




-Ey) 



{E-y - Ey) {E^^ - E^Y {Ej2 - Ey) 
-(-it)- 



{Ej - EyJ {E^.^ - {E^^ - Ey) 

^ \9l I 9l^ 9l '7772'?727'^77'' 



(B14) 



7i,'" ,75 



A?^ {nnc,ck) = 

1," 

= E 



E(-i)^-^ 



i=l 



^71 73 '?72 74 ^7375 '^71 74 ^71 7 ^7 



{E'j — -^71) {-^7 — E'y'Y {E^ — -^71) (-^7 — -^7') 



+ 



^-iEyt 



{Ej — -E'71 ) (-^71 ^ -^7') (^7 — E-y') (-E'71 ~ -^7') 



20" 



-i-it) 



{^Ej — Ejf^ (-^71 — -^7') 



(iz-y -^71) (-^7 -^7') 

2 



77 



„77l „7l7 
i/l i/l 



(B15) 
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AP (nnc, nk) = 

= E 



71 r-- 175 Li=l 



7ii72 



di{E[i,A]) 



Q — iE-yt 



n 



91 



'77173 ^77274 '^7375 ^771 74 '^717 '^757' 



(-E'7 ~ -^71) (-^7 ~ -^72) (-£'7 ~ -^7') 



+ 



(^^7 




- £"^2) (-£71 


-Eyf 












- E^^) {Ej^ 


- ) 


-Eyf 




e" 








- Ey) (S^i 


- Ey) {E^^ - 


-E,,f 




e" 


-iEyt 




(^7 


- iJy) {Ej^ 


- Eyf {E^2 


-Ey) 




e" 


-iEyt 





[E^ - Ey) {Ej^ - Ey) {E^2 - Ey) 

-i-it)- 



{E^ - Ey) {E^^ — Ey) {E^2 — Ey 



77i 7i7 

9i 9 1 ^772^7172^77' ' 



(B16) 



Finally, we calculate the Ap (nnn) by considering the two second decompositions, that is, its former three terms 



Ap {nnn,cc) = ^ 



7lj--- ,75 Li=l 



4 

n 



„7373 + l 

i/1 



Y\ '77fc7fc+2 j '^7174^7275^717^757' 



= E 



2Q-iEyt 



-\E.yt 



{E-y -E71) -£7') (-£7 ~ -£71) ~ -£"7') 



{E^ - E^^) (E^i - Ey) {E^ - Ey) {E^^ - Ey 



{E'y Eryty^ ("^71 ^'y' 



(£/-y -^71) (-^7 -^7') 



-(-it) 



g-iiJy* 



(£'-y Ej'^ (-£71 Ej') 



„77 „7 71 „7i7„77 



(B17) 



A4''' {nnn,cn) = 



7i>--- >75 



.1=1 



= E 



g— iE-,t 



J-J^7.7. + l 



Y\ ^7fc7fc+2 1 "^7174 ^7275^717*^757' 



Q—iEyt 



7y2 . (-^T ~ -^Ti ) (-^'^ ~ -^T^ ) (-^7 ~ Eyf {E^ - E^^ ) (i;^ - -E72 )^ (-£^7 ~ -^7' ) 



g— iS^t 



Q—iE^^ t 



{Ej — Ej^)^ {Ej - Ey^) {E^ - Ey) {Ej - E^j"^ {E^^ - E^^) {E^^ - Ey) 

■ + 



e~'^72* 



g-i£yt 



{E^ — -£72)^ (£^71 - -B72) (£72 - £7') (■E'7 ~ -^7')^ (-^71 ~ -^7') (-^72 ~ -£'7') 

^—iE~jt 



(£7 - £71) (£7 - £72) (-£^7 ~ -£'7') 



^77 ^7 7 ^7 '''(^J'''' rjj-^^^i rj^2^' . 



(B18) 
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AJ^ {nnn, nc) = 

7l.--- .• 



{Eh A]) 



Vlklk+2 1 '77l74'^7275'^7l7'^- 



757 



7i>72 
+ 



Q — iEy^ t 



{E^ - E^J {E^ — Ej^) {E^ — E^'f" {Ej - E^^) {E^^ — E^^) {E^^ — E^>)^ 



{E^ - Ejr) {E^^ - E^')"^ {Ery^ - Ey) {E^ - Ey)'^ {E^^ - Ey) {E^^ - Ey) 



^7 -^7 / V-^71 -^7 
-(-it) 



{Ery — Ej'^ {Eji — -^7') {Ej2 — E'y'^ 
while the fourth term has the third decomposition, that is 



„77 „7 71„7172 727 „ „ 

9l 9\ 9l fll '/77l'/772- 



Ap {nnn,nn,c) = 

7i>--- .■ 

= E 



di{E[j,4]) 



n 



„7j7j + i 

9i 



klk+2 J S.y^^5^^.y>6j^l 



7l>72,73 



Q—iE^t 



{E-y - -B^i) {E^ - E^^) {Ej - Ej^) 



[E, 


- E^^) [E^ - 


■ Ej^f {Ej - 




e~ 




{E^ 


- Ejif {E., 


- E^^) {E^ - E^^) 






lE-^^t 


(Ej 


- S^J^ (£^71 


~ ^72) (^71 ~ -^73) 




e~ 


i-^72* 


[E^ 


- E^^)^ {Ey^ 


- E^^) {Ej2 - E^^) 




e~ 





{E-y - -£73) (-E71 - £^73) (£72 - -£^73) 



(£7 ~ -^71) (E-y - £72) (-^7 ~ -^73) 

^ Si '''^ 9i^'^^9i^'^^ Si ^'''^772 ^71 73 ^77' • 



Ap {nnn,nn,n) = 

7ir-- >■ 

= E 



,i=l 



rfi(i^[7,4]) 



^7fc7i:+2 I '^7i7^757'%7' 



7ii72,73 



g— iS^t 



(-£7 ~ -^71 ) (-^7 ~ -^72 ) (Ej - -E73 ) (£'7 - £7' ) 
Q—iE^-^ t 



(E. 


-£,,)(£,, 




- £^.3) (£,,,j 


-£,0 


(Ej 


- £^2) (£'^1 


- -£^72) (-E'72 


- ^73) (-^72 


-£y) 


(^7 


~ -^73) (-^71 


- -^73) (-^72 


- -^73) {El3 


-£y) 


(^7 


- Ej') (£^71 


- £;y) (£'^2 


- Ey) (£;^3 


- Ey)_ 



+ 



^9'l^ 91 ^ 9l 91 ^ '7772 ^773 ^77''?7l73 ^7l7'^727' • 



(B19) 



(B20) 



(B21) 
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Now, all 15 contractions and/or anti-contractions in the fourth order approximation have been calculated out. 
In order to add the contributions from the fourth order approximation to the improved forms of lower order 

perturbed solutions, we first decompose AJ^ , which is a summation of all above terms, into the three parts according 
to e~'^^-*, (-ii)e"'^^-* and {-itye~'^-'^* /2, that is 



=A2''{e)+AJ''{te)+AJ''{fe). 
Secondly, we decompose its every term into three parts according to e~'^T*,e" 



e and e 



(B22) 



that is 



Ar'(e) = A2^'(e-'^-*)+Ar'(e-'^-i*)+Ar'(e-'^v*), 
Aj^'ite) = A2^'(te-'^-'*)+A2^'(te-'^-'i*)+A2^'(te-'^v*), 
Aj'^'it^e) = A2^'(t2e-'^-'*)+.l2'''(^'e-'^-i*)+A2^'(i2e-'^v*). 



(B23) 
(B24) 
(B25) 



Finally, we again decompose every term in above equations into the diagonal and non-diagonal parts about 7 and 7', 
that is 



le — ) = (e 
A^'it^ ~ 



D) 



-lE^.t. 



N), 



Al* (te-'^-'.*) = Al* (te-'^-.*;D) + {t 



-lE^.t. 



N), 



) = (t^e-'^^>*;D) + A2^ (t^e-'^^.*;N) 



(B26) 
(B27) 
(B28) 



where i?^. takes , and E-^i . 

If we do not concern the improved forms of perturbed solutions equal to or more than the fourth order one, we 
only need to write down the second and third terms in eq. ljB22p and calculate their diagonal and nondiagonal parts 
respectively. Based on the calculated results above, it is easy to obtain 



AV {te 



-iE-.,t 



-E 



9i 



111 |2 |yT">'2 |2 



[E^ - E^^y 

^71 72 



E 



771 1 2 |„772|2 



\9VT\9l 



72 



71,72 



71,72 ^-'^^ -^71) (-^7 ^72 y 
7l>72 



E 



[E^ ~ E-f^f {E^ - E^2 ) {E~f - E^^f' {E^ - E^^ ) 



E 



„77l„7l72 7273 737„ „ 

ill ill Hi ill '/772'/7l73 

[E^ - E^^) {Ej - E^^) {Ej - Ej^) 



(B29) 



Substituting the relation rjp-^fj^ 
of so-called revision energy G. 



3/31/32: 



using the technology of index exchanging and introducing the definitions 



= E 



71 



(B30) 



„77l „7l72 ^7273 „737„ 
.Vl t/l j/l t/l '/772 

{Ey - E^^ ) {Ey - Ej^ ) {E^ - E^^ ) 



E 



7li72,73 

we can simplify eq.(refA4gammaD) to the following concise form 



„77l „7l7„772 „727 

ffi 9i 9i ill 

71:72 (-^f ~ -^71) (^7 ~ -^72) 



(B31) 



(ie-'^-^D) = -(-it)e-'^-* 
Similar calculation and simplification lead to 

A]^' (te-'^-'*;N) = (-it)e-'^-* 



E 



^(2) 



71 (-^7 ~ -^71 ) 



/^(3) „77' 
<-^7 5l 



E 



G 



(2) 



(-E-y — Ej') (-^7 -^71) (-^7 ^ -^7') 



71 



(B32) 



(B33) 
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where 



„77l „7l72 ^727 

fi i/1 91 



[E^ ~ E^^) {E^ - E^^) 



(B34) 



For saving the space, the corresponding detail is omitted. In fact, it is not difficult, but it is necessary to be careful 
enough, specially in the cases of higher order approximations. 
In the same way, we can obtain: 



Ar'(fe-'^-*;D) = (-it)E,p p 

^1 ^71 



/-<(2) -iB.^ t 



(te-'^-i^N) - -Ht)E 



(^7 ~ -^71 ) (-^71 ~ ^7' ) 



771 7l7 

9{ 9i '777' 



(B35) 
(B36) 



(te-'^v*;N) = -(-i^e-'^v 



E 



(3) 



71 



{E^ - Ey) 



77 

9i 



E 



(2) 



g;i 



71 



(^/^ E.yf) (£'71 £7') 



771 7i7 

9i 9i VfY 



For i^e terms, only one term is nonzero, that is 



(-iGft)\-i 



2! 



since 



(B37) 



(B38) 



(B39) 

(B40) 
(B41) 



We can see that these terms can be merged with the lower order approximation to obtain the improved forms of 
perturbed solutions. 

2. 1=5 case 

Now let we consider the case of the fifth order approximation {I ~ 5). From ea. p2U|l it follows that the first 
decompositions of g-product have 2^ ~ 16 terms. They can be divided into 5 group 



(B42) 



where i indicates the number of 77 functions. Obviously 

A2'''{0;v)=Ar'{cccc), 



(B43) 



Al'' (1; ??) = AJ'^ (cccn) + AJ'^ (ccnc) + A}'' (cncc) + AJ'' (nccc), 



(B44) 



Al'' (2; 77) = AJ'' (ccnn) + AJ'^ (cncn) + AJ'' (cnnc) 

+Aj'^'{nccn) + AJ'^' (ncnc) + A^'^'(nncc), 



(B45) 
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AJ'^ (3; 77) = (cnnn) + (ncnn) + A'l'^ (nncn) + A'P (nnnc), 



(B46) 



AJ'' (4; 77) = AJ'' (nnnn). (B47) 

Here, we have used the notations stated in Sec. IVII 

By calculation, we obtain the A],'^ (0, 77) and every term of AJ'^ (1, rj) have only nontrivial first contractions and/or 

anti-contractions. But, we can find that every term of AJ^ (2, 77) can have one nontrivial second or third or fourth 
contraction or anti-contraction, that is 



(ccnn) = AJ"' 

(cncn) = AJ"' 

(cnnc) — AJ'^ 

(nccn) = AJ"^ 

(ncnc) — AJ"^ 



AT 
AT' 
AV 
AT 
AT 

Ap (nncc) = AJ"' {nncc, ckk) + AJ'^ {nncc,nkk). 
Similarly, every term of AJ'^ (3,??) can have two higher order contractions and/or anti-contraction: 



{ccnn, kkc) + AJ'^ (ccnn, kkn), 
{cncn, kc) -f Ap {cncn, kc), 
{cnnc, kck) + AJ'^ {cnnc, kck), 
{nccn, c) -I- Ap {nccn, n), 
{ncnc, ck) + A^ {ncnc, nk). 



(B48) 
(B49) 
(B50) 
(B51) 
(B52) 
(B53) 



AJ"^ {cnnn) 
Ap {ncnn) 
AJ"^ {nncn) 
AJ'^ {nnnc) 



AJ'^ {cnnn, kcc) -f- Ap {cnnn, ken) 
+Ap {cnnn, knc) -|- Ap {cnnn, knn), 
Ap {ncnn, kkc, ck) + Ap {ncnn, kkn, ck) 
+Ap {ncnn, kkc, nk) -t- Ap {ncnn, kkn, nk), 
A2j'^ {nncn, ckk, kc) + Ap {nncn, ckk, kn) 
+Ap {nncn, nkk, kc) + Ap {nncn, nkk, kn), 
AJ'^ {nnnc, cck) + AJ'^ {nnnc, cnk) 
+Ap {nnnc, nek) + Ap {nnnc, nnk). 



(B54) 
(B55) 
(B56) 
(B57) 



Moreover, their last terms, with two higher order anti-contractions, can have one nontrivial more higher contraction 
or anti-contraction: 

Ap {cnnn, knn) — Ap {cnnn, knn, kc) -\- Ap {cnnn, knn. kn), (B58) 
Ajp {ncnn, kkn, nk) — Ap {ncnn, kkn, nk,c) + A7P {ncnn, kkn, nk,n), (B59) 
Ap {nncn, nkk, kn) — Ap {nncn, nkk, kn,c) + Ap {nncn, nkk, kn,n), (B60) 

Ap {nnnc, nnk) — Ap {nnnc, nnk, ck) + Ap {nnnc, nnk, nk). (B61) 

In the case of Ap {nnnn), there are three second decompositions that result in 

Ap {nnnn) = Ap {nnnn, ccc) + Ap {nnnn, ccn) + Ap {nnnn, cnc) -f Ap {nnnn, ncc) 

-^Ap {nnnn, cnn) + Ap {nnnn, ncn) + Ap {nnnn, nnc) + AJ'^ {nnnn, nnn). (B62) 

In the above expression, from the fifth term to the seventh term have the third- or fourth- contraction and anti- 
contraction, the eighth term has two third contractions and anti-contractions: 



Ap {nnnn, cnn) = AJ'' {nnnn, cnn, kc) + Ap {nnnn, cnn, kn), 
Ap {nnnn, ncn) = Ap {nnnn, ncn, c) + Ap {nnnn, ncn, n), 
Ap {nnnn, nnc) — Ap {nnnn, nnc, ck) + Ap {nnnn, nnc, nk), 
Ap {nnnn, nnn) ~ Ap {nnnn, nnn, cc) + Ap {nnnn, nnn, cn) 

+Ap {nnnn, nnn, nc) + Ap {nnnn, nnn, nn) 



(B63) 
(B64) 
(B65) 

(B66) 
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In addition, AJ'* {nnnn, nnn, nn) consists of the fourth contraction and anti-contraction 

{nnnn, nnn, nn) — AJ'^ (nnnn, nnn, nn, c) + AJ'^ (nnnn, nnn, nn, n). 



(B67) 



According to above analysis, we obtain the contribution from the five order approximation made of 52 terms after 
finding out all of contractions and anti-contractions. 
Just like we have done in the I = 4 case, we decompose 



Ar {e) + Ar {te)+Ar (fe). 



(B68) 



where 



Ar (e) 



AJ'' (ie) = AJ^ (te-'^-'*) + AJ'' (fe-'^-i*) + AJ^ (te-'^-2*) + AJ'' (fe-'^v*), 
While, every term in above equations has its diagonal and non-diagonal parts about 7 and 7', that is 



(B69) 
(B70) 
(B71) 



^77 (e-ii^T.t) = ^77 (e-i^.,*;D) -f (. 



N), 



(B72) 
(B73) 
(B74) 



where E-y. takes Ery, E^-^, E^^ and E^i. 

If we do not concern the improved forms of perturbed solution more than the fourth order one, we only need to 
write down the second and third terms in ea. ljB68| l. In the following, we respectively calculate them term by term, 
and put the second and third terms in ea. (|B68|l together as A'P (te, i^e). 



{cccc; te, t^e) 



-(-it) 



3e 



3e 



iE^t 



{E^ - Ey) 
„-iE^,t 



(E^ - E^') 



2 {E.-E^r 



77 



77 



(B75) 



A^ {cccn; te, t^e) — ^ 



-(-it) 



2e-'^^ 



E^t 



{E^ — -^71 ) (-^7 — 



(-it) 



-(-it) 



(-it)2 



(Ey — i?^j)" {E^ — ^-f')"^ 



{E-y - E-y,f (E^, - Ey) 2! {E^ ~ E-y,) (E^ - E, 



7 



X 157''' I gr '7717' • 



(B76) 



A^ {ccnc; te, t^e) = ^ 



-it) 3 - (-it) 2 2 

{E^ — Ey-^ ) {Ej — Ej' ) (^7 — -^7' ) (-^71 ~ -^7' ) 



-(-it) 



2e-'^V* 



(-it)^ 



2! 



(i?-y i?^/ ) {E^^ Ejf ) 



,2 ,2 , 

X 5i'^ gF gr Viii- 



(B77) 



{cncc;te,t^e) = ^ 



-(-it) 



2e 



(i?-y -E-yj ) (-E-y — E-yl ) 



(-it) 



-lEyt 



{Ey — E^yi) (-^7 — -^7') 



(-it)2 



^71 / 

iE^t 



{Ey — Ej/) (-B71 ^7') ^' (-^7 ^ -^71) (-^7 ^ ^7' 



(B78) 
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AJ"' {nccc;te,t^e) = 



71 L 



-i-it) 



2e- 



3 - {-it) 3 

{Ej — -B-yi ) (-£'71 — -£'7') (-£'7 — E-y') ~ -£'7') 



-i-it) 



2! 



(^//y -^7') (-^71 -^7') 



^717 



51 ^771 • 



(B79) 



{ccnn, kkc; te, t^e) 



E 



7i>72 L 



-(-ii) 



-(-it) ^-"2 2 ~ (~^^) ~~3 

{Ef — {E^ — E^^) {Ey — E^J {E^ - E^^) 



+ 



(^7~-^7l) (^71^-^72) ^' i^l - E^l) (-£^7 --^72) 



x|5n'<?7''^5?'''5r^'577'- 



(B80) 



{ccnn,kkn;te,t'^e) = 



71,72 



(-it) 



{E^ - E^J {Ej - E^^) {Ej - E^>) 



+ (-«) 



(-£^7 ~ -^71)^ (-^71 ~ -^72) (■£'71 ~ -^7') 



X bi^'l gV^gT^^gT^ %i7'??772??77'- 



(B8I) 



{cncn,kc:tc,t'^c) 



E 



( 'i) 2c"'^^* ^ .^^ e~'^i'* 

(ii/7 — £"71) (-£'7 — Ej/)"^ {Ej — -£'71) (-£'7 — -£^7') 



+ (-ii) 



x\gT 



+ 



(-ii)^ 



{E^ — -^7') (-^71 — -^7') ^' (-^7 — -^71) (-^7 — -^7' 
2 



77 

gr 



77 

5i ^717' 



(B82) 



AJ"^ {cncn,kn;te,t^e) = 



-(-it) 



Q—iE^t 



7i>72 _ 

-(-it) 
-(-it)- 



{E^ — E^^) {Ej — E^^) {Ej - Ey) 
{E^ - E^J [E.f - E^i^f {E^ - Ey) 



+ 



{E-y - J^7i ) {Ej - E^^ ) {Ej - Ey ) 



2! {E^ — Ej^) {Ej — E^^) {Ej - Ey) 



I 5i ^7172 ^71 7' ^72 7' 



(BBS) 



Ag {cnnc, kck;te,t e) 



= E 



Q—iE^t ^—iE^-^t 

-it) 2 2 + 2 2 

(£/7 — -£^71) (-£'7 — -£'7') (-£^7 — -£'71) (-£'71 — -£'7') 

2 



+(-ii) 



{E.y — E^l) (-E71 -£'7') 



\gj 



111 I 



^717 



gV 



(B84) 



48 



{cnnc, knk;te,t^e) = 



7l>72 L 



i-it) 



{E^ - E^^ ) (Ej - Ej^ ) {Ej - Ey ) 



+{-it) 2- 

{E^ - Ej>) {E^^ - Ery,) - E^i) 



771 1 2 



^7.7 



Si ^7l7''?7l72'?772- 



(B85) 



Aj"' {nccn, c;te,t^e) 



E 



(-it) 



(-E'7 ~ ^71)^ (-^7 ~ -^72)^ ^ ' ^ (-^7 ~ -^71)^ (-^71 ~ -^72)^ 



(-^7 ^ ^72)^ (-^71 ~ -^72)^ 



1^7172 |2 ^771^7172^727^^^,. 



(B86) 



{nccn,n;te,t^e) = 



-(-it) 



(-E'7 ~ -^71) (-^71 ~ -^72) (-S71 ~ -^7') 



-(-it) 

(-^7 ~ -^72) i^-fi - -^72)^ (-£^72 - -E'7') 

^ l^i^'^^l 9i^^ 9i^^^ 9i^^ ^71 7' ^772 ^77'- 



(B87) 



Ag''' {ncnc,ck;te,t^e) 



= E 



71 



(-it) 



^ .^^ e-'^7'* 

(£/-y -^71) (-^7 -^7')^ (-^7 -^7') (-^71 -^7') 



-(-it) 



(-1/,)^ 



(i/zy — E^f^^ (-^71 — -^7') ^* (-^7 — -^7') (-^71 — -^7') 



7i7 

51 



77 

5i 



5i ^771 • 



(B88) 



(ncnc,nk;te,i^e) — 



7i>72 



-it) 



-(-it) 



-(-it) 



(-£^7 ~ -^7') (-^71 ~ -^7') (-^72 ~ -^7') 

(-£^7 - -E'7') (-E71 - -E^')^ (-B-y2 ~ -^7') 



(i;^ - Ey) {E^, - Ey) {E.,, - E.,,) 



(-it)2 



2! {Ej - Eji) {E^^ - Ey) (£72 ~ -^7') 
/ 2 ,2 , 
X g';^ g-^ 'ri"cyi'njj2Vji-y2- 



(B89) 



A5''' {nncc,ckk;te,t^e) 



= E 

7ii72 



-(-it) 



2e 



-iE-,t 



{Ej - EjJ {Ej - E^^) 



-lE-yt 



{Ej - E^^f {E^ - Ej^f 



(-it)^ 



+(-it) 3 ■ 91 2 

{E^ - EyJ {Ej^ - Ej^) 2. {E^ - E^J {Ej - Ej^) 



(B90) 



49 



{nncc,nkk;te,t^e) = 



i-it) 



71)72 L 



(-^7 ~ -^72) (-^71 ~ -^72) (-S72 - ^1') 



(-E'7 ~ -^7') (-^71 ~ -^7') (-^72 ~ -^7') 



^ Si^''' g'T'^ gl^'^^ gl^'^ ^77i7''?772 ^77' 



(B91) 



^5''' {cnnn, kcc; te, t^e) 



E 

7i>72 



-(-it) 



-(-«) 



{Ej - Ej^) {E^ - Ej^) 



(-£^7 ~ ^lif (-^7 ~ -^72)^ 



2 



2! 



(-^7~-^7i) i^l - ^12) 



x\grrgrgr'grs,^'- 



(B92) 



A5''' {cnnn, ken; te,t^e) = ^ 



(-it) 



{Ej - EjJ {E^ ~ E^^) {Ej - Ejr) 

(^7 - Eflf (^71 ~ -^72) (-^71 ~ -^7') 



^|^77.|^^772^727l^7.7^^^^,^^^,. 



(B93) 



[cnnn, knc; te, t^e) = ^ 

7l>72,73 
-(-it) 

-(-it)- 



-(-it) 



(-E'7 - -2^71 ) (-^7 " -^72 ) {^1 - -E'73 ) 

^—'lEyt 



{Ej - Ej^ ) {E^ - Ey^ ) {E^ - Ej^ ) 



+ 



{E-Y - -Et,i) (-E^ - E^^) {E^ - E^^) 
{-itf e-'^^* 



2! (-E'7 — -£'71) (-£"7 ^ -^72) (-^7 ^ -^7 



^ I fll '''^ S'l^^'771 72 '?7l 73*^77' ■ 



(B94) 



{cnnn, knn, kc; te, t^e) = 



71.72 



(-it) 



{Ej - E^^ ) (£^ - E^^ ) (E^ - Ey 
{Ej - Ej')^ {Ej^ - Ey) (£^2 - £y ) 



/ 2 , 

X gV ??772%i7' 



(B95) 



aI ' {cnnn, knn, kn; te, t e) 



7i>72,73 



e ' Iffl'*' I gT^ gj gj ^773 ^77'^7l72^7l73^7l7'^^727' 

{E^ - E^J {E^ - E^^) {E^ - E^^) {E^ - E^) 



(B96) 



50 
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{nncn, nkk, kc; te, t^e) = 



7ii72 



-(-it) 



{E^ - E^J {E^^ - E^J {E^^ - Ey) 



+(-ii) 



{E^ - Ej') {E^^ - Ey)^ {E^2 - Ey) 



(B104) 



Ap {nncn, nkk, kn, c; te, t^e) = 

->72' 

-(-it) 



7i>7273 L 



(-it) 



Q—iE^t 



{E^ - EjJ{E^- E^^ ) {E^- E^^) 



{E^ - E^^f {E^^ - E^^) {E^^ - E^^) 



(B105) 



AJ'^ {nncn, nkk, kn, n; te, t^e) 

^ I 'y273 |2 771 7172 727' 

_ ^ -^-^ S ^2 I ^772%73'^77''?7l73^7l7'^737' 

^.■r:^^, (-^7 " -^72) (-E'71 ~ -^72) (-^72 - -^73) (-^72 - -^7') 



(B106) 



Ag (nnnc, cck; te, t e) 



E 



71)72 L 



-(-it) 



(E"^ - -E-yi)^ (E^ - E^a)^ (E^ - E^-i)^ (-Bt. - 



-(-it) 



^—'lEy^ t 



+ 



{E^ - EjJ {E^^ - E^^) 2! _ E^^) {E^ - E^^) 



X 1^771 |2 ^771 ^7172^727^^^,_ 



(B107) 



{nnnc,cnk;te,t^e) = 

71)72,73 

-(-it) 
-(-it)- 



-(-it) 



[E-y - E^J {E^ - Ej^) {E^ - Ej^) 

Q—iEyt 



{Ey - E^J [E^ - E^i^f {E^ - E^^) 



+ 



{E-Y - -E^i ) {Ej - -E^2 ) {E-Y - -E^a ) 

(-it)2 



Q—iEyt 



2! {E^ - E^^) {Ej - E^^) {Ej - Ej^) 

X \9l'^^\ ^'''^fl'l^''' %i 73*77273 ^77'- 



(B108) 



Ap {nnnc, nek; te, t^e) = 



7i>72 L 



-(-it) 



+(-ii) 



{Ej - E^J {Ej^ - Ey^) {E^^ - Ey) 

(Ef - -By) {E^^ — Eyf' {Ej^ - E^i) 
2 / 

„77i„7i72 727 „ „ , 
91 9l 9l '7772 '/77- 



(B109) 



52 



{nnnc,nnk, ck;te,t^e) = 



7i>72 



i-it) 



{E^ - E^^) {Ey - E^^) {Ej - Eji) 



2- 

{Ef - -E7') (£'71 - E^>) {Ej2 - Ej') 



(BllO) 



AJ'^ {nnnc, nnk, nk; te, t^e) 



^737 



^77 gJ^g'J'^ ^772'7773^77'^7i73^7i7'^7273 



(-2^7 ~ -^7') (-^71 ~ -^7') (-^72 ~ E-r') (-^73 ~ -^7 



(Bill) 



A2'' {nnnn,ccc;te,t^e) 



{-it) 2 2 + (~^^) 2 2 

(i/^ — -^71) (^7 — -^7') (-^7 — -^71) (-^71 — -^7') 



(i/^ — -^7') (-^71 — -^7') 



/\ 2 



1^771 |2 (^^771^717- j ^77' 



(B112) 



Ag'*^ {nnnn,ccn;te,t'^e) = ^ 



(-ii) 



71,72 L 



{Ej-EjJ {E^ - E^^) {Ej - Ey) 



Q—iE^^ t 

{E^ - E^^f {E^^ - £^7 J (-5^71 ~ -^7') 



X(57^^5r^^ffr^)57^^5?"^77'^727'- 



(B113) 



Ag''' {nnnn,cnc;te,t^e) = ^ 



i-it) 



7l>72 L 



(£7 - £71) (£7 - £72) (£7 - £7') 

(£7 - £7') (£71 - £7') (£72 - E7/) 



XlflZ-'^ff?^"') (57^^57="') 57"'^7i72- 



(B114) 



vlg''' {nnnn,ncc;te,t^e) = 



-(-it) 



(£7 - £7 J (£7^ - £"72) (£7^ - £7/) 



(£7 - £7') (£71 - £7') (£72 - £7') 

) 



(B115) 



AJ'^ {nnnn,cnn,kc;te,t'^e) = 



(-it) 



(£7 - £71 ) (£7 - £72 ) (£7 - £7' ) 
(£7 - £7')^ (£71 - £7') (£72 - £7') 



x{9t'''9r) {9r9r'')9r'v,^,.- 



(B116) 
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A'^'^ (nnnn, cnn, kn; te, t^e) 



-iB-t „77l „7l72 _727 „773 „737 



• 5l 51 .9l 51 .9l '?77'^7i73^7273^7l7''7727' 



(-£^7 ~ -^71) (-^7 ~ -^72) (-£'7 ~ -^73) (-^7 ~ -^7') 



(B117) 



Ag'*' {nnnn,ncn,c;te,t^e) = ^ 

71)72 73 

+(-ii) 



(-it) 



(-£^7 ~ ^71)^ (^7 ~ -^^72) (-^7 ~ -^73) 



Q iEy^ t 



(^7 ~ -^71) (^71 ~ -^72) (^71 ~ ^73) 
^ Iffl'' I 5l ' 9l ^ 9l ^ '?772'7773'^77'' 



(B118) 



Ag''' {nnnn, ncn, n; te, t^e) 



i-it) 



C lE^^t gjll g'ril2 g-y273 gjall gjlf 7] ^ ^ , 7] ^ ^ ^ 7] ^ ^ ^ fj ^ ^ I 7] ^ ^ J ' 



(B119) 



AJ'^ {nnnn,nnc, ck;te,t^e) = 



(-it) 



{E^ - ) {E^ - E^2 ) (-^7 - -^7' ) 



+(-it) ^ 

{E^ - Ej') {E^^ - E^f) {E^^ - Eji) 



(B120) 



Ap {nnnn, nnc, nk; te, t^e) 

„-i-E^/t„7'72 7273 737' 771 „7l7' ,n r, n n 
_ (it) '/77 '/772'/773'/7l72'/7l73 



(B121) 



^Ig"*^ {nnnn, nnn, cc; te, t^e) = 



7i>72 



(-it) 



Q-iEyt 



{Ej - E^^ ) {Ej - E^^ ) (£^ - Ey 

e-'^y* " 

{E~f - E~f>) {E^^ - Eryt) {Ej2 - E^,) 



xgr92''9r'9r9rV,-y.V,.Y 



(B122) 



A'l'^ {nnnn, nnn, cn; te, t^e) 

^ [f^ ^ ^1 ^1 '''' 9l''''' 9l^ ^772^7l73^7l7'^727'^737' 



7li72,73 



(-2^7 ~ -^71) (-^7 ~ -^72) (-^7 ~ -^73) (-^7 ~ -^7') 



(B123) 



A'l'^ {nnnn, nnn, nc; te, t^e) 

(-^7 - ^7') (-2^71 - ^7') (-^72 - Ej') (^73 - ^7') 



(B124) 



{nnnn, nnn, nn, c; tc, t^e) 

e ' ^ gj^ gj^gj^ 9i''' 9i ''^772 ^773*771 73*771 74 '77274 '^77' 



7i>72,73,74 



(-2^7 ~ -^71 ) (-^7 ~ -^72 ) (-^7 - ^73 ) (-^7 ~ -^74 ) 



(B125) 



54 



AJ'^ {nnnn,nnn,nn,n;te,t'^e) = 



(B126) 



Based on above 52 contraction- and anti contraction- expressions, we can, via the reorganization and summation, 
obtain 



where 



71)72 



71,72,73,74 



{E^ - Ey^) (Ey - E^^) 



{E^ - EjJ {E^ - E^^) {Ej - E^^) {Ej - E^J 



{Ej - E^^ ) {Ej - E^^ ) 

^771 ^7172^727 J^^, + {-iG(^h)6^,, 



- E 



„77l „772 „7l 7 „72 73 „737 

i/l Hi Hi Hi Hi 



+ 



„77l „772 „7l7 „72 73 „73 7 
iJl ill i/l i/l i/l 



+ - 



{Ej - Ej^)'^ {Ej - Ej^) {Ej - Ej^) {Ej - EjJ {E^ - E^^)'^ {Ej - Ej^) 

„77l „772 ^717 „7273 „737 
i/l i/l i/l i/l i/l 



(-E'7 ~ -^71) (-^7 ~ -^72) (-£'7 ~ -^73) 



(B127) 



(B128) 



71 (^7 ~ -^71 ) 

+(-iGf i) ^ 



rgr'grs-y 



71,72 (-^7 ~ -^71) (-^71 ~ -^72) 



^771^7172^727^^ 



(B129) 



e~'^72* 



71,72 (-^-Y ~ -^72) (-E^7i ~ -^72 
A5(te-'^v*,D) = 



.^771^7172^727 J 



77 



(B130) 
(B131) 



— iE t "71 

l5(ie-'^^*,N) = (-iGWi)^-f- + (-iGft)^ 

^77/ 7j 



e-ii;,tg77lg7l7'^77' 
(-2^7 ~ -^71) (-^7 ~ -^7') 



-(-iG(^^i)E 

71 

+(-iG?)i) E 



e-iE,t^77l^7l7^77 



+ 



e-iB,t^77l^7l7^77 



e-'^^'g7"^g7^^^gr"'%72r?77' 

(-E'7 ~ -E'71) (-^7 ~ -^72) (-£^7 ~ -^7') 



(B132) 



A5(ie-'^-i*,N) = -(-ii)E 

71 



gCTe-ii;-,itg77lg7l7'^^^, 
(-^7 ~" -^71) (-^71 ~ -^7') 

G(?c"'^^i*.g7-^\g7^-^-.g7--^'7j7,yr;,y 

(-^7 ~ -^71) (-^71 ~ -^72) (-E'71 ~" -^7') 



(B133) 



71,72 



g(2),-ii.-„t^^77lg7l72g727'^^^^^^^, 

{E^ - E^^) {E^^ - E^^) {E^2 - ^7') 



(B134) 
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[E^ - Ey ) 

71 



71 

3-iS,,*^7'7l^7l7'^77' 



e-ii;ytg77lg7l7'^77' 

{E^ Ejl) (-E'yj — Eyi) 



e-^E-r'^gJ ''^gf'' g-p 



[Ey — E~f') {E,y-^ — ^7') (-^7 — ^7') (-^71 ~ -^7') 

e-'^-'*57"^57^^^5r"'^7i7'^77' 



{Ey Ey' ) {Ey^ — EyI ) {Ey2 " -E'y 



(B135) 



For the parts with t'^e, we have 



2-iE-.t 



2! ^""7 '-'7 "77' 



A5(i'e-'^-i*,D) = A5(t"e-'^-2*,D) = ^(i'e-'^v^, D) = 



J.2 _ — iE-y^ t 



(B136) 
(B137) 



2! ^ 

A5(i'e-'^-2*,N) = A5(i2e-'^v*,N) = 0, 



yl5(t'e-'^v*,N) = - 



2! 



2 e-'^v* 



77 



(B138) 
(B139) 
(B140) 



It is clear that above diagonal and nondiagonal part about AJ^ (te) and AJ'^ (te) indeed has the expected forms 
and can be merged reasonably to the lower order terms in order to obtain the improved forms of perturbed solutions. 

3. 1 = 6 case 

Now let we consider the case of the sixth order approximation (/ = 6). From ea. p2(J|l it follows that the first 
decompositions of ^-product have 2^ = 32 terms. Like the I = 5 case, they can be divided into 6 group 



1=0 



where i indicates the number of rj functions. Obviously 

A'/{0;v)^Al'''{ccccc), 



(B141) 



(B142) 



(i;'7) 



AJ^ (ccccn) + A^'^ (cccnc) + AJ'^ (ccncc) 
+A2'^ (cnccc) + Ap (ncccc), 



(B143) 



Ar (2;r/) 



Ap (cccn) + A'P (ccncn) + AJ^ (cnccn) + A^ (ncccn) 
+A'f (ccnnc) + A'P {cncnc) + AJ"' (nccnc) + A'P (cnncc) 
+A'P {nance) + A'l'^ {nnccc), 



(B144) 



Ar (3;^) 



A'P (ccnnn) + AJ^ [cncnn) + A'l'^ (cnncn) + AJ^ (cnnnc) 
+A1'^ (nccnn) + AJ'^ (ncncn) + Aq'^ (ncnnc) + AJ'^ (nnccn) 
+^g^ (nncnc) + AJ^ (nnncc), 

Aq'^ (4; rj) — AJ^ (cnnnn) + A^"' (ncnnn) + AJ^ (nncnn) 
+^g^ (nnncn) + Af {nnnnc) 



(B145) 



(B146) 
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AT (5;??) = ^T (nnnnn). 
Furthermore considering the high order contraction or anti-contraction, we have 



(B147) 



AT {ccnnn 
AT {cncnn 
AT {cnncn 
AT {cnnnc 

I 

AT {nccnn 
AT {ncncn 
AT {ncnnc 
AJ"' (nnccn 
AT {nncnc 
AT {nnncc 



AT (cccnn) 


= AT 


(cccnn, kkkc) + AT {cccnn, kkkn), 


(B148) 


AT (ccncn) 


= AT 


(ccncn, kkc) + AT {ccncn, kkn), 


(B149) 


AT (ccnnc) 


= AT 


(ccnnc, kkck) + AT (ccnnc, kknk), 


(B150) 


AT (cnccn) 


= AT 


{cnccn, kc) + AT {cnccn, kn). 


(B151) 


AT {cncnc) 


= AT 


{cncnc, kck) + AT {cncnc, knk). 


(B152) 


AT (cnncc) 


= AT 


{cnncc, kckk) + AT {cnncc, knkk). 


(B153) 


AT {ncccn) 


= AT 


{ncccn, c) + AT {ncccn, n), 


(B154) 


AT {nccnc) 


= AT 


{nccnc, ck) + AT {nccnc, nk). 


(B155) 


AT (ncncc) 


= AT 


{ncncc, ckk) + AT {ncncc, nkk). 


(B156) 


AT (nnccc) 


= AT 


{nnccc, ckkk) + AT {nnccc, nkkk) 


(B157) 



= AT {ccnnn, kkcc) + AT {ccnnn, kkcn) + AT {ccnnn, kknc) 

+AT {ccnnn, kknn, kkc) + AT {ccnnn, kknn, kkn), 
= AT {cncnn, kkkc, kck) + AT {cncnn, kkkc, knk) + AT {cncnn, kkkn, kck) 

+AT {cncnn, kkkn, knk, kc) + AT {cncnn, kkkn, knk, kn), 
= AT {cnncn, kckk, kkc) + AT {cnncn, kckk, kkn) + AT {cnncn, knkk, kkc) 

+A2'^ {cnncn, knkk, kkn, kc) + AJ'^ {cnncn, knkk, kkn, kn), 

= AT {cnnnc, keck) + AT {cnnnc, kcnk) + AT {cnnnc, knck) 

+A2'^ {cnnnc, knnk, kck) + AT {cnnnc, knnk, knk), 
/ / / 

= AT {nccnn, kkkc, ck) + AT {nccnn, kkkc, nk) + AT {nccnn, kkkn, ck) 

+AT {nccnn, kkkn, nk, c) + AT {nccnn, kkkn, nk, n), 
= A2^ {ncncn, ckc) + AT {ncncn, ckn) + A^'^ {ncncn, nkc) 

+AT {ncncn, nkn, c) + {ncncn, nkn, n), 
= AT {ncnnc, kkck, ckk) + AT {ncnnc, kkck, nkk) + AT {ncnnc, kknk, ckk) 

+AT {ncnnc, kknk, nkk, ck) + AT {ncnnc, kknk, nkk, nk), 
= AT {nnccn, cc) + AT {nnccn, cn) + Aq ' {nnccn, nc) 

+AT {nnccn, nn, c) + AT {nnccn, nn, n), 
= AT {nncnc, ckkk, kck) + AT {nncnc, ckkk, knk) + AT {nncnc, nkkk, kck) 

+AT {nncnc, nkkk, knk, ck) + AT {nncnc, nkkk, knk, nk), 

III 

= AT {nnncc, cckk) + AT {nnncc, cnkk) + AT {nnncc, nckk) 
+Aq"' {nnncc, nnkk, c) + AT {nnncc, nnkk, n) 



II I I 

AT {cnnnn) = AT {cnnnn, kccc) + AT {cnnnn, keen) + AJ'^ {cnnnn, kcnc) 

+AT {cnnnn, knee) + AT {cnnnn, kenn) + Aq^ {cnnnn, kncn) 

I I 

+A2'^ {cnnnn, knnc) + AJ^ {cnnnn, knnn). 



(B158) 
(B159) 
(B160) 
(B161) 
(B162) 
(B163) 
(B164) 
(B165) 
(B166) 
(B167) 



(B168) 



Aq'^ {cnnnn,kcnn) 
AJ'^ {cnnnn,kncn) 
{cnnnn,knnc) 



A'g' {cnnnn, kcnn, kkc) + A'^'^ {cnnnn, kcnn, kkn), 
AJ'^ {cnnnn, kncn, kc) + AJ^ {cnnnn, kncn, kn), 
A'P {cnnnn, knnc, kck) + Af {cnnnn, knnc, knk). 



A^"^ {cnnnn, knnn) 



Aq^ {ncnnn) 



= Ag'*' {cnnnn, knnn, kcc) + A'^'^ {cnnnn, knnn, ken) 
+Ag^ {cnnnn, knnn, knc) + {cnnnn, knnn, knn, kc) 
^-A'P {cnnnn, knnn, knn, kn). 

Aq'^ {ncnnn, kkcc, ckk) + AJ'^ {ncnnn, kkcc, nkk) 
+Aq"' {ncnnn, kkcn, ckk) + Ap {ncnnn, kknc, ckk) 
+AJ'^ {ncnnn, kkcn, nkk) + AJ'^ {ncnnn, kknc, nkk) 
+Aq'^ {ncnnn, kknn, ckk) + Aq' {ncnnn, kknn, nkk). 



Ap {ncnnn, kkcn, nkk) 
Ap {ncnnn, kknc, nkk) ■- 
Ap {ncnnn, kknn, ckk) - 

Ap {ncnnn, kknn, nkk) 



= A^'^ {ncnnn, kkcn, nkk, c) + A^ {ncnnn, kkcn, nkk, n), 
= Ap {ncnnn, kknc, nkk, ck) + AJ^ {ncnnn, kknc, nkk, nk), 
— Ap {ncnnn, kknn, ckc) + A^ {ncnnn, kknn, ckn), 

Ap {ncnnn, kknn, nkc, ck) + A^ {ncnnn, kknn, nkc, nk) 
+A2'^ {ncnnn, kknn, nkn, ck) + AJ'^ {ncnnn, kknn, nkn, nk, c) 
+Aq^ {ncnnn, kknn, nkn, nk, n) 



Ap {nncnn) = A^ {nncnn, ckkc, kck) + A^ {nncnn, ckkc, knk) 

+Ag^ {nncnn, ckkn, kck) + Af {nncnn, nkkc, kck) 
+AJ'' {nncnn, ckkn, knk) + Ap {nncnn, nkkc, knk) 

{nncnn, nkkn, knk). 



+j4g {nncnn, nkkn, kck) 



■AP 



AJ'^ {nncnn, ckkn, knk) 
Ap {nncnn, nkkc, knk) 
Ap {nncnn, nkkn, kck) 



Ap {nncnn, ckkn, knk, kc) + A^ {nncnn, ckkn, knk, kn), 

I I 

Ap {nncnn, nkkc, knk, ck) + Aq^ {nncnn, nkkc, knk, nk), 
Ap {nncnn, nkkn, kck, c) + AJ"' {nncnn, nkkn, kck, n). 



AJ^' {nncnn, nkkn, knk) 



Ap {nncnn, nkkn, knk, cc) + Ag ' {nncnn, nkkn, knk, cn) 
+Ap {nncnn, nkkn, knk, nc) + A^ {nncnn, nkkn, knk, nn, c) 
+Ap {nncnn, nkkn, knk, nn, n) 



Ap {nnncn) 



Aq ' {nnncn. cnkk, kkn) 
AJ'^ {nnncn, nckk, kkn) 
Ap {nnncn, nnkk, kkc) 



AJ"^ {nnncn, cckk, kkc) + AJ"^ {nnncn, cckk, kkn) 
+Aq'^ {nnncn, cnkk, kkc) + A^ {nnncn, nckk, kkc) 
+Ap {nnncn, cnkk, kkn) + Ap {nnncn, nckk, kkn) 
-\-Ap {nnncn, nnkk, kkc) + Ap {nnncn, nnkk, kkn), 

v4g ' {nnncn, cnkk, kkn, kc) + AP {nnncn, cnkk, kkn, kn), 
Ag^ {nnncn, nckk, kkn, c) + Ap {nnncn, nckk, kkn, n). 



Ap {nnncn, nnkk, ckc) 



■ Ap {nnncn, nnkk, nkc), 
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{nnncn, nnkk, kkn) = A^'^ {nnncn, nnkk, ckn, kc) + A^^ {nnncn, nnkk, ckn, kn) 

+Ag^ {nnncn, nnkk, nkn, kc) + Aq'' {nnncn, nnkk, nkn, kn, c) 
{nnncn, nnkk, nkn, kn, n). 



(B187) 



Aq' {nnnnc) = Al'^ {nnnnc, ccck) + AJ^ {nnnnc, ccnk) + Ag'*' {nnnnc, cnck) 

r / t 

+Ap {nnnnc, neck) + {nnnnc, cnnk) + AJ^ {nnnnc, ncnk) 

{nnnnc, nnnk), 



+Aq {nnnnc, nnck) 



■AT 



(B188) 



A^ {nnnnc, cnnk) 
Ap {nnnnc, ncnk) 
Ag'*' {nnnnc, nnck) 



AJ"' {nnnnc, cnnk, kck) + AJ^ {nnnnc, cnnk, knk), 
Ap {nnnnc, ncnk, ck) + AJ"' {nnnnc, ncnk, nk), 
Ap {nnnnc, nnck, c) + AJ'^ {nnnnc, nnck, n). 



(B189) 
(B190) 
(B191) 



Ap {nnnnc, nnnk) 



AP {nnnnc, nnnk, cck) 



-AP 



g ^/t/t/t/to, jvjvjvrv, uoA, J -p -^g {nnnnc, nnnk, cnk) 
+Ap {nnnnc, nnnk, nek) + AJ'^ {nnnnc, nnnk, nnk, ck) 
+Aq^ {nnnnc, nnnk, nnk, nk). 



(B192) 



Ap {nnnnn) = A^"^ {nnnnn, cccc) + AJ"^ {nnnnn, cccn) + AJ^ {nnnnn, ccnc) 



{nnnnn, cncc) + AJ'^ {nnnnn, nccc) + AJ'^ {nnnnn, ccnn) 
{nnnnn, cncn) + A^ {nnnnn, cnnc) + A^ {nnnnn, nccn) 
{nnnnn, ncnc) + A^ {nnnnn, nncc) + AJ'^ {nnnnn, cnnn) 
{nnnnn, ncnn) + A^ {nnnnn, nncn) + A^ {nnnnn, nnnc) 



+Aq' {nnnnn, nnnn). 



(B193) 



Ap {nnnnn, ccnn 
Ap {nnnnn, cncn 
Ap {nnnnn, cnnc 
Ap {nnnnn , nccn 
Ap {nnnnn, ncnc 
Ap {nnnnn, nncc 



AJ'^ {nnnnn, ccnn, kkc) + Ap {nnnnn, ccnn, kkn), 
Ap {nnnnn, cncn, kc) + Ap {nnnnn, cncn, kn), 
Ap {nnnnn, cnnc, kck) + Ap {nnnnn, cnnc, knk), 
Ap {nnnnn, nccn, c) + Ap {nnnnn, nccn, n), 
AJ'^ {nnnnn, ncnc, ck) + Ap {nnnnn, ncnc, nk), 
Ap {nnnnn, nncc, ckk) + Ap {nnnnn, nncc, nkk). 



(B194) 
(B195) 
(B196) 
(B197) 
(B198) 
(B199) 
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Aq'* {nnnnn, cnnn) = AJ'^ {nnnnn, cnnn, kcc) + Ag'"' {nnnnn, cnnn, ken) 

+Aq'^ (nnnnn, cnnn, knc) + AJ^ (nnnnn, cnnn, knn, kc) 

+Aq'^ (nnnnn, cnnn, knn, kn), 
Ap (nnnnn, ncnn) = AJ'^ (nnnnn, ncnn, kkc, ck) + AJ'^ (nnnnn, ncnn, kkc, nk) 

+A'P (nnnnn, ncnn, kkn, ck) + A'^ (nnnnn, ncnn, kkn, nk, c) 

+A'P (nnnnn, ncnn, kkn, nk, n), 
(nnnnn, nncn) — (nnnnn, nncn, ckk, kc) + AJ^ (nnnnn, nncn, ckk, kn) 

+Ag'^ (nnnnn, nncn, nkk, kc) + Ag'^ (nnnnn, nncn, nkk, kn, c) 

+Aq''' (nnnnn, nncn, nkk, kn, n), 
AJ"^ (nnnnn, nnnc) = A'P (nnnnn, nnnc, cck) + AJ"^ (nnnnn, nnnc, cnk) 

+Aq'^ (nnnnn, nnnc, nek) + AJ^ (nnnnn, nnnc, nnk, ck) 



A^fP (nnnnn, nnnn) 



+Ap (nnnnn, nnnc, nnk, nk), 

Ap (nnnnn, nnnn, ccc) + A^ (nnnnn, nnnn, ccn) 
+A2'^ (nnnnn, nnnn, cnc) + A^ (nnnnn, nnnn, ncc) 
+A2'^ (nnnnn, nnnn, cnn) + A^ (nnnnn, nnnn, ncn) 
+Ap (nnnnn, nnnn, nnc) + A^ (nnnnn, nnnn, nnn). 



(B200) 



(B201) 



(B202) 



(B203) 



(B204) 



(nnnnn, nnnn, cnn) — A^ (nnnnn, nnnn, cnn, kc) + A^ (nnnnn, nnnn, cnn, kn) , 
Ap (nnnnn, nnnn, ncn) — AJ^ (nnnnn, nnnn, ncn, c) + A^ (nnnnn, nnnn, ncn, n) , 
Ap (nnnnn, nnnn, nnc) — Ap (nnnnn, nnnn, nnc, ck) + Ap (nnnnn, nnnn, nnc, nk) , 



(B205) 
(B206) 
(B207) 



Ag^ (nnnnn, nnnn, nnn) = AJ'^ (nnnnn, nnnn, nnn, cc) + AJ^ (nnnnn, nnnn, nnn, cn) 

-\-Ap (nnnnn, nnnn, nnn, nc) + A^ (nnnnn, nnnn, nnn, nn, c) 

+Ap (nnnnn, nnnn, nnn, nn,n). (B208) 

Thus, we obtain the contribution from the six order approximation made of 203 terms after finding out all of 
contractions and anti-contractions. 

Just like we have done in the / = 4 or 5 cases, we decompose 



AJ'' = AP (e) + AJ-' (te) + A^ (fe) + AJ'' (t^e) 
= Ap'(e,te)+Ap'(t\t'e), 



(B209) 
(B210) 



To our purpose, we only calculate the second term Ap (t'^e, t^e) in ea. (|B210p . Without loss of generality, we decompose 
it into 



AP (t\t^e) = AP (t\-'^-'\fe-'^-'*) + AP (t^e-'^-'^\t^e-'^-'^') 



+Ar (fe 



While, every term in above equations has its diagonal and non-diagonal parts about 7 and 7', that is 

AP' (t^e-'^-'') = Ap' (th-'''-n';B) + Ap' (fe-'^-'^';N), 
Ap' (t'^e-'^-'^') = Ap' (t^e-'''-"';B) + Ap' (t^e-'''-'^';'N). 



(B211) 

(B212) 
(B213) 



where Ej. takes E^,E^^ and E^i. In the following, we respectively calculate them term by term, and we only write 
down the non-zero expressions for saving space. 
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AQ{ccncn,kkn;t^e,t^e) = ^ 



(-i^)^ c-'^^i* 1.97^^1^ \gr'\'grgr'v-n.V,yV^.Y 
2 (S^ - E^^) {E^^ - E^^) {E^^ - Ey) 

(_i^)2 c-i^-'* 



71,72 



AQ{ccnnc,kkck;t^e,t^e) = ^ 



77 



„77l „7l7 



Ae{cnccn,kc;t^e,t^e) = ^ 



71 2 {Ej - E^j"^ {Ej - Eyf 



77 



771 7i7 

gigi 



2 (-^7 — -^71) (-^7 — -^7') 



(B223) 



(B224) 



(B225) 



{-■^tf e-' ^^* Ig?-^^!^ |g?"1^g7'^^9^"'^77'^7.7.^7 l 
(-E'7 ~ -^71 ) (-^7 ~ -^72 ) (^7 ~ -^7') 



(B226) 



AQ{cncnc, kck;t'^e,t^e) = 



3! (jEJ^-E^J {E^ - Ej^) 



2 - E^^) {E^ - E'^J 



-\E.,t 



Ae{cncnc, knk;t^e,t'^e) = 



X 157^'!'' IS'l'^'l^»?7i72'^77'- 



3! {Ej - E^^) {Ej - Ej^) {E^ - Ej^) 

Q — iE-yt 



2 (i;,- 


- E^^) {E^ - E^^) {E^ - E^^f 






2 (E,- 


- E^^) {Ej - E^^f {Ej - E^^) 






2 (E,- 


- E~f^f [E^ - E^^) [E^. - Ey^) 


1^771 |2 1^772 


Idl^ 1 ^7172 ^7173 ^7273 '^77'" 



A^ {cnncc, kckk; t^e, t^e) = ■ 



gr 



77 

gV 



g^gr 



i^Ej — E'yiy^ (-^71 — -^7')'^ 



{-it) 

1 

A6(ncccn, c; t^e, t^e) = - ^ ^ 



71,72 



(ncccn, n; t^e, f^e) = — 



(-it)^C-'^->i*|.97-^M'|gr"1%772V 

2 {E^ - Ej^f {Ej^ - E^^y 

(_i^)2 e-ii^.1* |gr^1%7"^gr"'^77'^772 ^727' 

2 {E^ - E^^) {E^^ - E^^f' [E^i - Ey) 



AQ{nccnc,ck;t^e,t^e) = 



i-it) 



2 e ^ 



77 



^717 



„77l „7l7 



71 



2 {E^ — E^' ) {-^71 — ) 



(B227) 



(B228) 
(B229) 

(B230) 

(B231) 

(B232) 
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AQ{nccnc,nk;t^e,t^e) = 



i-it) 



2 e-'^-v'* 



7l7 

9i 



9l 9 1 ^77''?772^7i72 



(B233) 



AQ{ncncc,ckk;t^e,t^e) = 



7i>72 L 



Q — iEyt 



{E^ - Ej^ ) {E^ - Ej^ ) 



(-ii)^ 



1^77. |4 1^7172 |2^^^^^^^, 



(B234) 



AQ{ncncc,nkk;t^e,t^e) = — 



9r 



{E^ - EjJ {Ej^ - E^^) {E^^ - Ey) 



(B235) 



AQ{nnccc,ckkk;t^e,t^e) = 



71 



(-ii)^ 



+ 



2 (-£'7 — -^7') (-^71 — -^7') 



77 



^771^717 



(B236) 



2„ +3„\ _ \ " 



AQ{nnccc,nkkk;t e,t e) = 



(-it) 



2 e~'^^'* 



„72 7 



^77 gj 7 rjj^irjjj^rjj-^j^ 



-Y, .-7, (-^7 ~ ^7') (-^71 ~ -^7') (-^72 - ^7') 



(B237) 



Ae{ccnnn,kkcc;t^e,t^e) = 



2 (-S-y — -£^71) (-^7 — E^')"^ 



2 (£^7 — Ej')"^ (-^71 — -^7') 



57" 



„77 „7 71 „7l7„77 

9i 9i 9i 9i 



(B238) 



AQ{ccnnn,kkcn;t e,t e) 



Er-i/^S -i_E^t I 77l|2 77l„7l72 727„77' ,„ , 
( 11) e ^ ig;^ I g;^ ff^ ^7l7'^727' 

71 ,72 ^ (-^7 - ^71 ) (-^7 ~ ^72 ) (-^7 ~ ^7' ) 



(B239) 



Ae{ccnnn,kknc;t e,t e) 



i-it) 



2 e-'^v* 



77 



^77 gj'^gj^gj'^ V'y'ri '^772 



7^0^; ^ (-^'>' " -^7')^ (-^71 ~ -^7') (-^72 - ^7') 



(B240) 



Aeiccnnn,kknn,kkc;t'e,t^e) = 



(-ii)' e-'^^'|g7^^|^g7^^g7^^-g7-^--gf^ry,^,7?-,,^3J, 



77 



2 (-£^7 — -^71 ) (-^7 — -^72 ) (-^7 — 



13) 



(B241) 



A(, icncnn, kkkc, kck; t^e, t^e) = 



77 

57 



^771^717 



71 



2 (i;^ - E^^f iE^ - Ey f 



(B242) 



Ae{cncnn,kkkc,knk;t^e,t^e) = 



(_i^)2 e-'^^*|g7^^|^ 



77 

57 



51" 57 " ^7172^717' 



71,72 



(B243) 



63 



71,72 



{E~f - E^^ ) {Ej - E^2 ) (-^7 ~ ) 



(B244) 



iE^t I „77l 1 2 



{cncnn, kkkn, knk, kc; t^e, t^e) = 



(_it)2e \g{ 



772 727 „ „ , 
fl al '/7i72'/7i7 



2 - i;^^ ) (i;^ - E^^ ) {Ery - Ej> ) 



(B245) 



AQ{cncnn, kkkn, knk, kn; t^e, t^e) 



E 

71 ,72 ,73 



(_ii)2 e-ii^.* 1^771 1^1^77. 

I ffl ffl ^77'^7l 72^71 73 ^7l7'^7273 ^727' 
(-£^7 ~ -^71) (-^7 ~ -^72) (-E'7 ~ -^73) (-£'7 ~ -^7') 



(B246) 



Ae{cnncn,kckk,kkc;t^e,t^e) = 



7i>72 



2 - E^if {Ej - E^^ ) 



77 



(B247) 



AQ{cnncn,knkk,kkc;t^e,t^e) = 



(_i^)2 ^-iE^t |g77i|2 1^772 |2 [g^^^a |2 ^^^^^^^^^^^^^^^^ 

2 {Ej - Ej^ ) (i;^ - Ej^ ) - Ej^ ) 



(B248) 



AQ{cnnnc,kcck;t^e,t'^e) = 



71 



(-ii)^ 



2 (-£"7 — -^71) (-^7 — -^7')^ 



+ 



2 (-£/7 — -^7') (-^71 — -^7') 



77 



771 7i7 

9i 9i 



(B249) 



A(,{cnnnc,kcnk;t e,t e) 



2„ i3„\ _ 



(-it) 



2 e~'-^7'* 



77 



„727 



^77 ^7 7 Tj^^^Tj^^^^ 



-71 ,-72 ^ (-^7 ~ -^7')^ (-^71 ~ -^7') (^72 - ^7') 



(B250) 



AQ{cnnnc,knck;t^e,t^e) = 



77 



772 727 „ „ , 

i/l fl '/7i72'/7i7 



7i>72 



2 - E^^ ) (iJ-y - £'-^2 ) {E^2 - Ej 



\2 ■ 



(B251) 



A^{nccnn,kkkc,ck;t^e,t^e) = 



(-it)2e-'^-'* 



^717 



„77 „7 7i„7i7„77 

9i 9i 9i 9i 



71 



2 (-E'7 — -E'7') (-^71 — E-y') 



(B252) 



Ae{nccnn,kkkc,nk;t^e,t^e) = 



7i7 

5i 



^7 7 ^7 7 ^7 7 jy^^^jy^^^ 



71,72 



(-^7 ~ -^7') (-^71 ~ -^7') (-^72 - -^7') 



(B253) 



(_ii)2 e-'^-*|57"N' 



^^717 



aVgr ^77' 



(■E'7 ~ -^71)^ (-^71 ~ Ey)'^ 



(B254) 



7ii72 



(£^7 ~ -^71) (-^71 ~ -^72) (-^71 ~ Ey) 



(B255) 
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Ae{ncncn,nkc;t^e,t^e) = — 



{-it) 



2 e '-^71* 



1^7172 I grJ''''igrT'iT rij^ir]jj2Tlj'y''n'y2j' 



{E^ - EjJ {E^^ - Ej^) (E^i - Ey) 



(B256) 



{-itf e-'^^i* \gr f \9r' f igr^ i' %72^773%273^7 



Aq (ncncn, nkn, c;t e, t e) 



^ 2 

71.72,73 



{E^ - E^J {Ej^ - Ej^) {E^^ - Ej^) 



(B257) 



AQ{ncncn, nkn, n; t e,t e) 



■ E 

7l>72,73 



I 3l 5l ^772'7773 ^77''77273'7727'^737' 
(£'7 ~ -^71) (-^71 ~ -^72) (-^71 ~ -^73) (-^71 ~ -^7') 



(B258) 



AQ{ncnnc,kkck,ckk;t^e,t^e) = 



77 



(B259) 



Ae{ncnnc,kknk,ckk;t^e,t^e) = 



71,72,73 



e-ii^.* IgJ-'^f |g7^3|2 |g7^^^|%^^^^^,^^,3^^^^^j^^, 

2 - E^-,) {E^ - E^^) {E^ - E^^) 



(B260) 



AQ{nnccn,cc;t'^e,t^e) = 



71 



(-it)2e-'^7* |g77l|2^77g7'7lg7l7g77' 
2 (-E'7 ~ -E'71) (-E'7 — -^7') 



(B261) 



A6(nnccn, cn; t^e, f^e) = 



(-ii)^ e-'^-* |57"l'ff7^^flr"^5r"97'''^7i7' 



77727' 



71 ,72 ^ {Ej - EjJ {E^ - E^^f {E^ - Ey) 



(B262) 



n2 e^'-^^'* 



57" 



^717 



^771^717 



Ae(nncnc,ckkk,kck:t'^e,t^e) = ^ ^ ^ ^ 

^ ^ 2 {E,-EyfiE,,-Eyf 



(B263) 



Ae{nncnc,ckkk,knk;t^e,t^e) = 



(-i<) 



2 e"'^7'* 



ff7" 



^727 



771 71 7 
9l 9l '7772'?7l72 



V, ,7, 2 (-B7 -Eyf {E^^ - Ey) (£^2 - Ey) 



(B264) 



AQ{nncnc,nkkk,kck;t e,t e) 



2„ +3„\ _ 



(-ii) 



2e '^7'* 



^717 



„72 7 



g,77 gJ 7 77^^2'?7i 72*777' 



71,72 ^ (-^7 - ^7') (-^71 ~ ^I'f (-^72 - -2^7') 



(B265) 



AQ{nncnc,nkkk,knk,ck;t e,t e) 



2p f3^x ^ V- 



i-it) 



2 e-'^v* 



77 



gr 



71,72 ^ (-^T - Eyf {E^^ - Ey) (£^2 - Ey) 



(B266) 



AQ[nncnc, nkkk, knk, nk; t^e, t^e 



(-it) 



2 '^V* 



^727 



^737 



^77 gi 7 ri'yj2'ni^3'nji''n^i^2'n^i^3'n^2^3 



^ 2 

71,72,73 



(-2^7 - -E'7') (-^71 ~ Ey) (£'72 - £'7') (-B73 - £7') 



(B267) 



65 



AQ{nnncc,cckk;t^e,t^e) = 



71 



2 (-£'7 — -^71 ) (^7 — ^-y')"^ 



{Ej — Ejf^ (-^71 — -^7') 



77 



„77 „7 71 „7i7„77 



(B268) 



AQ{nnncc,cnkk]t e,t e) 



i-it) 



2 e 



5?" 



„77l „7l72 72 7„77 „ ,„ , 
91 9l 9l 9l '/7i7''/727' 



^ 2 

71.72 



{E^ - E-yJ (i;^ - E^^) {E^ - Eyf 



(B269) 



Ae{nnncc,nckk;t^e,t^e) = 



2 c 



-iE-,l 



Hi 



Hi Hi ' Wi 'Hi' ih-.i'h 



{Ej - Ey) {Ej^ - Ey) {E^^ - Ey) 



(B270) 



AQ{nnncc,nnkk,c;t'^e,t^e) = 



{-it? e-^^* \9r'f9r9r'9r'9rv,-y.v-y.,sS. 



77 



(-^7 ~ -^71) (-^7 ~ -^72) (-^7 ~ -^73) 



(B271) 



Ae{cnnnn,kccc;t e,t e) = 



(-i^)2 e-i^.* \9r'f9r'9j'''9r9r' 

^ 2 {E,-E,,)\E,-E,,f 



(B272) 



Ae{cnnnn,kccn;t^e,t^e) = 



7i>72 



e-ii^.* \gj'>^ f 9r'9r'9r9r'v,.YV,.Y ^ 

(-£^7 ~ -^71)^ (-^7 ^ -^72) (-^7 ~ -^7') 



(B273) 



/ :j-\2 1^771 |2 77' 7'72 727 „77' ,„ 

Ae{cnnnn,kncc;t\t'e) = ^ " '■'^^ ' ^"^^^^^^ 



(B274) 



AQ{cnnnn,kcnn,kkc;t'^e,t^e) = 



77 



7ii72i73 



2 (E^ - i;^! ) {E^ - Ej^ ) {E^ - £'^,3 ) 



(B275) 



Ae {cnnnn, kncn, kc; t^e, t^e) = ^ 



-it) 



2 e"'-^^* 



77 



771 7172 727 77 , 



7l>72 



2 (E^ - ) {E^ - ) (£^ - £7' )^ 



(B276) 



Ae{cnnnn, kncn, kn; t e, t e 



E 

71 >72 ,73 



73 7 



2 (£7 - £71) (£7 - £72) (£7 - £73) (£7 - £7') 



(B277) 



Ae{cnnnn,knnn,kcc;t e,t e) 



i-itf 0-'^-.* |57"^|'57"=yr^^5r^^5r^^773^7i72'57r 



^ 2 

71,72,73 



(£7 - £71 ) (£7 - £72 ) (£7 - £73 ) 



(B278) 



AQ{cnnnn, knnn, knc: t^c, t^c) 

( — e 1,97'' I 9l^ 9l '^^ gj^^ gj ^?7773T?7i72?]7i73?]7i74'772 74'^77' 



E 



71,72,73,74 



(Ef - £71) (£7 - £72) (£7 - £73) (Ej - £74) 



(B279) 



66 



Ae{ncnnn,kkcc,ckk;t^e,t^e) = 



i-it) 



2 e~'-^7'* 



77 



7i7 

9i 



2 (-^7 — -^7') (-^71 — -^7') 



(B280) 



Ae{ncnnn,kkcc,nkk;t^e,t^e) = ^ 



(-ii) 



2 e^^^T'* 



Si'' 5l ffl '* 3l '?77i^772 



(B281) 



{ncnnn, kknc, ckk; t^e, t^e) = 



i-it) 



2 e-'^y* 



57" 



^717 



^77 gj 7 rjj^-^TJj-^-y^ 



71,72 ^ (-^7 -^I'f (^71 " -^7') (-^72 - -^7') 



(B282) 



Ae{ncnnn,kknc,nkk,ck;t^e,t^e) = ^ 



(-if) 



2 e-i^y* 



71 7 

5i 



77 7 72 727 77 „ „ 
i/l i/l i/l i/l '/77l'/7l72 



(B283) 



Ae{ncnnn, kknc, nkk, nk: t^c, t^e) 

2 



= E 



i-it) 



2 e"'^^'* 



7i7 

5i 



.97''' 9i'^9i'^9i'^ '7771 ^772 '7773 ^7172 ^7173 



2 (£^7 - ) (i?7i - iJ-y/) {Ej^ - Ey) [E^^ - E^>) 



(B284) 



A6{nncnn,ckkc,kck;t^e,t^e) = ^ 



(-it)^e~^^^* |g7"M'lg7"l'|g7^"l'V ^ 



(B285) 



{nncnn, ckkc, knk; t^e, t^e) = 



(_it)2 e-ii^.* |g7^=^|^g7^^g7^^^gr''^5f "'^7173^7 



71 ,72 ,73 ^ (-E'7 - -^71 ) (-^7 - -^72 ) ^ (£^7 - -^73 ) 



(B286) 



AQ{nnncn,cckk,kkc;t^e,t^e) = 



'iE^t |„77i 



1.97' 



.9i' 



771 717 
9i 9i 



2 (i;^ _ Ej^f {E^ - Ej') 



(B287) 



AQ{nnncn,cckk,kkn;t^e,t^e) = 



71,72 



(_i^)2 ^-iE.,t |ff7^^|^97^^g7^'^^97-^ff7-^;,^,,r77,7, 

2 (i;^ - E^j"^ {Ej - £^72 ) (£7 - £7/ ) 



(B288) 



iE^t I „772 |2 



AQ{nnncn,cnkk,kkc;t e,t e) = 



(_if)2e '^^*|g7 



57" 



„77l„7l7„ „ 
fl al '/7172 '/727 



2 (E.^ - ) (E.^ - Ej^ ) (£7 - Eji ) 



(B289) 



(-_7^|2 ^iB^t I 772|2 77' 7'7i 717 77' 

Ae{nnncn,cnkk,kkn,kc;t\t'e) _ ^ ^ l-9i I 5i 5i 9i 9i '^7172^727' 



71,72 ^ (-E'7 - -^71) (^7 - -^72) (-^7 - ^7')' 



(B290) 



AQ{nnncn, cnkk, kkn, kn; t e, t e 



E 

71 ,72 ,73 



( it)^ e |ff7'^^l 57'*' ff7 '*' g7 ''^57''' ''77l73''77l7''^7273''7727''?737' 

2 {E.y - Ej^ ) (iJ^ - E^^ ) - £"^3 ) - £.^/ ) 



(B291) 



67 



AQ{nnnnc,ccck;t^e,t^e) = 



i-it) 



2 e~'^^^'^ 



7i7 

9i 



„77 „! 7i„7i7„77 

yi yi i/i i/i 



2 (-^7 — -E'-y' ) (-^71 — -£/-y' ) 



(B292) 



Ae{nnnnc,ccnk;t^e,t^e) = 



77 7 71 71 7 77 

9i 9i 9i 9i '7772'?7i72 



7li72 



(B293) 



Ae{nnnnc,ncck;t^e,t^e) = 



i-it) 



2 e"'^^'* 



gT^ gl^ gl'' gl ^771 '7772 



{E-y - £"7') (£'71 - -E7') (£72 - £7') 



(B294) 



AQ{nnnnc,ncnk,ck;t e,t e) 



2„ +3„\ _ 



i-it) 



2 e^'-^-r'* 



gr 



gV gr^gT'^gr r,-y,.v,,. 



71,72 ^ (-^7 ~ ^7') (^71 - -^7') (-^72 - ^7') 



(B295) 



Ae{nnnnc,ncnk,nk;t e,t e 



2 e-'^v* 



^737 



gj77 gj'^gj'^g'j'^ ^771 '?772 ^773 ^7173^77273 



71,72,73 



(£7 - £7') (£71 ~ -^7') (-^72 ~ -£'7') (-^73 ~ -^7') 



(B296) 



{nnnnc, nnck, c; t^e, t^e) 



e-i^.* |g7^- 1^ g7^\97^"^gr"^g?^'77i73^77' 



1 (-£^7 ~ ^71 ) (-^7 - -^72 ) ^ (-2^7 ~ ^73 J 



(B297) 



^ f-lV^S -iS^i I 77l|2 771 7172 7273 737 J 

^6(nnnnc,nnnfc,ccfc;t2e,i3e) = ^ '^^ ' ^72^7173^ 

71,72,73 (-B7--B71) (£7 - -B7J (£7 - -E73) 



(B298) 



AQ{nnnnc, nnnk, cnk; t^e, t^e) 

(—if) e ' 1'* l^i''^ I gj^ gj^gj^ gj '*^^772'?773 ^7i73^7i74 ^7274 ^7374*^77' 



E -2 

71,72,73,74 



(-2^7 - -E'71 ) (-E'7 - £72 ) (-£^7 ~ -^73 ) (-^7 ~ -^74 ) 



(B299) 



AQ{nnnnn,cccc;t'^e,t^e) = 



(_i^)2e-i^.*(g7-rig7i72g727)2 ^^ 

2 {Ej-E^J {Ej-E^^) 



"77 



(B300) 



^ C-if^2 p-i£;,t„77i„7i72 727 „77i„7i73 737„ X , 

Ae{nnnnn,ccnc;t\t^e) = ^ ( if) e ^ g, g, g, g, g, v^^^.d^j' _ ^gg^^^ 

71,72,73 (-2^7 ~ -^71) (-^7 ~ -^72) (-E^7 ~ ^73) 



^ f_;/-^|2 -iE^t 77l„7l72 727 773 7372 727 I 

A6(nnnnn,cncc;f2e,f3e) = ^ ^^^^ ''^^ ^\ ^7173^7' _ (^302) 

71,72,73 ^ (-£^7 ~ -^71) (-^7 ~ ^72) (^7-^73) 



_ C_:^\2 -i_E^t 771 7172 727 773 7371 717 , 

A^{nnnnn,cnnc,kck;t\t^e) = ^ ^7273')77' _ (^303) 

71,72,73 ^ (-^7 ~ -^71) (-^7 ~ ^72) (-E'7 ~ -^73) 



68 



Ae{nnnnn, cnnc, knk; t^e, t^e) 

\ '■'') ^ Ul Ul Hi Hi .Vl iJl '/7l73'/7l74'/72 73'/72 74 "77 (B304) 

71 ,72 ,73 ,74 ^ (^7 - -^71 ) (^7 " -5^72 ) (-^7 " "^73 ) ("^7 " "^74 ) 

Based on above non-zero 91 terms and vanishing 112 terms with t^e^f'^e faetor parts in ah of 203 contraction- and 
anti contraction- expressions, we can, via the reorganization and summation, obtain the following concise forms: 

-^-^y ^gV'gr^^i'- (B305) 

(2)\^„-iE,,t 



A6(t^e-'^-*,D) = ^-J^Y^h-Er^-ET^aV^ar^ii'. (B306) 



AQ{t^e-'^^'\V)) = 0. (B307) 



71 



12 ICg'fo-'^^.' 



Their forms arc indeed the same as expected and can be merged reasonably to the lower order terms in order to 
obtain the improved forms of perturbed solutions. 



2 
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